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Abstract. Local bifurcation theory typically deals with the response of a degenerate but isolated 
equilibrium state or periodic orbit of a dynamical system to perturbations controlled by one or 
more independent parameters, and characteristically uses tools from singularity theory. There 

arc many situations, however, in which the equilibrium state or periodic orbit is not isolated 
but belongs to a manifold S of such states, typically as a result of continuous symmetries in the 
problem. In this case the bifurcation analysis requires a combination of local and global methods, 
and is most tractable in the case of normal nondegeneracy, that is when the degeneracy is only 
along S itself and the system is nondegenerate in directions normal to S. 

In this paper we consider the consequences of relaxing normal nondegeneracy, which can 
generically occur within 1-parameter families of such systems. We pay particular attention to 
the simplest but important case where dim S = 1 and where the normal degeneracy occurs with 
corank 1. Our main focus is on uniform degeneracy along S, although we also consider aspects 
of the branching structure for solutions when the degeneracy varies at different places on S. The 
tools are those of singularity theory adapted to global topology of S, which allow us to explain 
the bifurcation geometry in natural way. In particular, we extend and give a clear geometric 
setting for earlier anals^tical results of Hale and Taboas. 
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1. The general problem 

1.1. Introduction 

Let Fq : M TM be a smootb(|] vector field on a smooth (possibly infinite-dimensional) manifold 
M, with the property that the flow generated by Fq has a smooth d-dimensional manifold S of 
equilibria with finite c? > 0. In other words, vanishes everywhere on S. We now perturb _Fo to 

F,:^F^ + eT+0{\e\'') (1.1) 

where T is 'a smooth vector field on a neighbourhood of S in M and where £ is a real parameter, 
and we ask what can be said about zeros of F^ close to 5' in M for e sufficiently small but nonzero. 
In particular, we ask which points G S* are such that a branch of equilibria emanates from them 
as e varies away from zero. 

More generally we may consider a multiparameter deformation Fg, where £ S M"^, as in j2^[57] 



or [TTl [T5] for example. In this case the term eT in (1.1 1 is interpreted as J2i=i^i-^i ^'^^ ^ 
collection {Ti , . . . , !Fq} of smooth vector fields, and 0(|eP) means terms of degree 2 or higher 
in the components of e. Here it is useful to think in terms of equilibrium branches as e moves 
away from the origin in a given direction in parameter space R'^, and we discuss this more carefully 
below. Of course, it is not clear a priori that branching behaviour depends only on J- and not 
on the 0(|ep) terms, although it is a key aspect of our results that often under suitable generic[|] 
assumptions on !F this is indeed the case. 

In many situations where bifurcation from a manifold arises in applications the manifold S 
is normally non-degenerate with respect to _Fo i meaning that at each point x Cz S the tangent 
space TxS (which automatically lies in kerT^Fo) in fact coincides with herT^Fo, so that T^Fq 
is as non-degenerate as possible under the given geometric assumptions. Bifurcation under these 
conditions has been well studied: see [2 , [T3] , [TB] , [TO] , [IH] , [37] • A standard technique is to show by 
a global Lyapunov- Schmidt reduction along 5* that for small e there is a smooth manifold close 
to S (viewed as the image of a small section of the normal bundle NS of S in Af ) such that the 
'normal' component of F^ vanishes on S^. This reduces the problem to a study of the 'tangential' 
component of F^, which corresponds to a tangent vector field F^ on S via the projection NS — > S. 
The sought-for zeros of F^ close to S then correspond to the zeros of the vector field F^ on S itself|jj] 

In the absence of further information the number of zeros of -F^ can often be given a lower 
bound by topological means such as index theory: see the techniques of [THlUnj for example. In the 
variational case whe re F^ is the gradient field of a real-valued function the same may be assumed for 



(see Section 2.1 1 and the usual tools are Lyusternik-Schnirelmann category theory in general or 
Morse theory if critical points are assumed to be nondegenerate: see for example . [5] . [55] . . [T5] 
and [35], [3] II] for applications to Hamiltonian systems. With further information on symmetries 
[H] > [31] , [7] , [37] or on the actual perturbing field T more specific results can be obtained as in 
[10|[IIl[ni[13| and [ll[15lll6]. 

I By smooth we mean C'°° . However, it is clear that many of our results hold with only a few orders of 
differentiability. 

§ In our context a generic condition or assumption is one which holds for an open dense set of the maps or objects 
under consideration. Making this precise requires formal machinery of transversality in jet spaces and techicalities 
which we mostly prefer to avoid: in our cases the interpretation of genericity is straightforward. 

II For a more careful description see Section [2] 
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The purpose of this paper is to investigate the consequences of relaxing the normal non- 
degeneracy condition by supposing that T^S is proper subspace of kerT^^Fo , that is 

dim kerT^Fo = k + d (1.2) 

for all X Cz S, where A; > 1 does not depend on x . We call this normal degeneracy with constant 
corank k. Our most detailed results are for the simplest (but important) case when fc = 1, but some 
of the methods and observations are applicable also to the general case. 



1.1.1. The group orbit case. There is a natural situation in which manifolds of equilibria arise 
automatically, namely where the manifold S is the orbit of the action of a Lie group T under 
which Fq is equivariant. In particular this occurs when Fq is the gradient of a F-invariant function 
fo : M M.. In problems formulated on function spaces of periodic functions, there is a natural 
circular symmetry defined by the phase-shift, and in mechanical systems natural symmetries arise 
from coordinate invariance. Translational symmetry can also play a key role in bifurcation theory 
of differential equations on the real line: see [30 for example. Problems that we consider, in which 
the perturbed system Fg^ does not share the symmetry ( F-equivariance) of Fq , go under the general 
heading of forced symmetry-breaking. 

For a typical group orbit we expect normal non-degeneracy, but if Fq is a member of a smooth 1- 
parameter family of vector fields Fq each of which is F-equivariant then, since (for square matrices) 
dropping rank by 1 is a codimension-1 phenomenon, we would expect to find isolated values of the 
parameter 7 G M for which the corresponding equilibrium manifold is normally degenerate with 
constant corank 1 . The constancy of corank in this case is a consequence of the equivariance. 

1.1.2. Non-constant corank. In contrast with the above, non-degeneracy may fail in that 
dimkerT^j-fo > dimS* only at isolated points x €z S . This does not arise in the group orbit 
case, since the group action implies that the behaviour at all points of S is 'the same'. Analysis of 
this type of degeneracy requires different techniques, which we touch upon in Section |5] See also 

m- 

1.1.3. Sources of the problem. The break-up of families of periodic orbits of Hamiltonian systems 
is a classical problem, first treated in a modern spirit in the fundamental papers of Weinstein |38j 
and Moser [32 . Abstract treatment of bifurcation from a manifold was pioneered by Hale [2B] , with 
applications to planar forced ordinary differential equations, and extended by Hale and Taboas [55] 
to the case of normal degeneracy. Further important techniques and results relating particularly to 
symmetries were set out by Dancer |19L 1201 1^ and Vanderbauwhede |37j . again supposing normal 
nondegeneracy. Applications of bifurcation from a manifold were exploited by Ambrosetti et al. [21 14] 
in the context of Hamiltonian systems and also more generally [5]. Meanwhile, a general geometric 
framework for multiparameter bifurcation from a normally nondegenerate manifold was formulated 
by the first author [TnidllllS], and customised tools adapted to families of periodic orbits of ordinary 
differential equations were systematically refined in a series of papers by Chicone [HI [TSJ [T2] • It is 
the work of Hale et al. and Chicone that is the main inspiration for the current paper. 

1.2. Outline of the approach 

An overview of our approach in this paper is as follows. First we apply Lyapunov-Schmidt reduction 
to enable us to focus on the core of the problem. As a consequence of normal degeneracy, the 
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reduced problem still involves a fc-dimensional normal component of the vector field, as well as the 
automatic c?-dimensional tangential component. We would like to 'unfold' this degenerate {k + d)- 
dimensional vector field, using e as a multi-parameter, but we are unable to use standard tools of 
unfolding (deformation) theory directly as the degeneracy along S causes the local problem to have 
infinite codimension. Therefore we make appropriate generic assumptions about the behaviour of 
the deformation term J- along the manifold S. With this, we can in principle describe generic 
geometry of such deformations, at least for small k . In the fundamental case k = 1 we are able to 
give a rather complete picture. 

In Section|2]we formulate the initial reduction process, and in Section[3]we show how to adapt 
the standard notion of versal deformation to our context, giving particular attention to the simplest 
nontrivial low-dimensional case. In Section [4] we combine global geometry with local algebra to 
describe generic branching structure, and in Section [5] we study local branching structure in closer 
detail. Finally, in Section |6] we discuss three areas of application of the methods set out in this 
paper. 



2. Lyapunov- Schmidt reduction 

Let NS be a normal bundle for S in M, that is a vector sub-bundle of TgM := TM\s orthogonal 
to TS (using a Riemannian structure i' on M) . A neighbourhood of the zero section of NS can be 
taken to parametrise a tubular neighbourhood of 5 in M via the exponential map associated to 
and so then we may without loss of generality regard the whole bifurcation problem as taking place 
on a neighbourhood W of the zero section of the bundle NS. 

Let H he a, connection on NS, that is a (smooth) sub-bundle of T NS orthogonal to (vectors 
tangent to) the fibres of NS: thus 

T^NS = e T^N, ^H^®N, 

for all w lying in the fibre N^ of NS over x £ S. This enables us to define a linear bundle map 
TT : T NS NS as projection onto the second factor in each tangent subspace T^^NS. 

Given a vector field F : W C NS TNS the composition nF : W ^ NS satisfies 
TrF{W n Nx) C Nx for each x £ S. Thus ttF is a 'vertical' vector field on W obtained by 
projecting F parallel to the 'horizontal' subspaces determined hy H. If it; e n N^ we write 
DF{w) G Lin{Nx, Nx) to denote the 'vertical' derivative of F (i.e. the derivative of F\Nx) at w. 

Let F^ : M ^ TM , e e M« be a vector field as in Section [l] with Fo{x) = for all x € S. 
As above, we replace M by C NS. Since T^S C keiT^Fo for all x £ S, the constant corank 



hypothesis (1.2 1 implies that at the origin 0^ of N^ the kernel := ker DFq(Ox) is a /c-dimensional 



linear subspace of N^ for every x G S* . We now make the following assumptions: 



(Al) The subspaces {K^} together define a smooth /c-dimensional vector sub-bundle K of NS. 
(A2) There is a smooth codimension-fc vector sub-bundle L — {Lx} of NS complementary to K in 

NS ; thus Nx^ Kx® Lx for all x E S. 
(A3) There is a smooth vector bundle decomposition NS = P ® R such that Rx is the range of 

DFq{Qx) and dimPx = k for all x £ S. 
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Remark. In finite dimensions the smoothness of K holds automatically, given the constant 
dimension k, since the kernel of a matrix is the orthogonal complement of the range of its transpose: 
thus if DFq{Ox) varies smoothly with x then a basis for K" orthogonal to the rows of DFo{Ox) can 
also be chosen to vary smoothly. The construction of L and P is straightforward. 

Smoothness of K also holds in many naturally-arising infinite-dimensional cases where DFo{Ox) 
is a differential operator, so that its kernel consists of the solutions ^{t) to a /c-dimensional system 
of linear differential equations; construction of L and P is also often easy in such cases. 

If DFo{Ox) is self-adjoint (as for a gradient vector field, for example) then we may take = Rx 
and Px = Kx- 

Let 5*0 denote the (image of the) zero section of the bundle K. Of course 5*0 = S, but the 5*0 
notation will focus attention on 5 as a submanifold of K . Applying Lyapunov-Schmidt reduction 
globally along a neighbourhood of 5*0 in K we obtain the following result. 

Proposition 2.1 There exists a neighbourhood U of Sq in K with the property that for all 
sufficiently small e there is a unique smooth bundle map : U ^ L such that 

(i) cr"|5o = and 

(ii) the R-component of F^ vanishes precisely at the points lying on the graph C NS of . 

Proof. Let r denote the bundle projection NS = P ® i? — s- R, and consider the composition 
r o ttF^ : W ^ R. When e = the vertical derivative of this smooth map at 0^; G Nx (identified 
with X (z So) is by definition of i? a linear surjection with finite-dimensional kernel Kx. Hence by 
the Implicit Function Theorem (IFT), given xq G Sq there is a neighbourhood V of xq in K and 
£o > such that for all |e| < Sq and x ^ V O Sq the solution set to r o irF^ in the fibre Nx is 
the intersection with Nx of the graph of a unique smooth bundle map ay : V ^ L . See Figure 
[T] The construction of on a global neighbourhood U of Sq from local sections ay then follows 
by uniqueness of ay on V and the local compactness of S. □ 

Since the (c? -I- fc)-dimensional manifold is precisely where the vertical component of the 
vector field F^ has zero i?-component, it follows that the zeros of F^ are the zeros of the vector field 
on If^ obtained by projecting F^ to TS(BP. (Note that this projected vector field is not necessarily 
tangent to .) The problem therefore distils to that of solving the reduced problem 



where 1 denotes the identity map on TM = TNS. Here if w = {u,a'^{u)) G Nx we identify 
with TxS through the natural bundle projection TNS TS, so that 1 — r o tt is projection on 
TS®P. 

Remark. If F is of class then F is also for 1 < p < co, w. 
2.1. The gradient case 

In the important special case when F^ is the gradient V/e of a real- valued function (with respect 
to the Riemannian structure i^) then the zeros of correspond to the zeros of a gradient field 
on U. To see this, take R = L and P — K (see earlier remark) and observe that \i u & U and 



Fe{u) := (1 



r ott) F,{u,a%u)) ^OeTS®P, 



ueU C K 



(2.1) 



fe{u) ■= fe{u,a^{u)) then 
V/e(w) 



(l-ro7r)V/e(M,cr^(u)) 



(2.2) 
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= graph{ay] 



Figure 1. The geometry of Lyapunov-Schmidt reduction at xq £ 5o • Here A' is the normal 
bundle to 5 = So in M and Kj,^ = kei DF{xo). 

there being no contribution to V/e(w) from the second argument of because r o ttV/^ vanishes 
identically on C/^ by construction. 



2.2. Interpretation of the reduced problem 

The reduced vector field Fe ■ U C K ^ TS ® P satisfies 

Fo{0^) = for all a; e 5* (2.3) 

and 

DFoiOa:) = for all x e S, (2.4) 
since Fq ~ and cr" = on 5*0, and also (1 — r o tt)DFo{Ox) = by definition of r. We can therefore 
regard Fg, as arising from a family of smooth singular map germs 

parametrised (globally) by cc € 5 which is then perturbed (locally) by e e M'' for |e| small. After 
shrinking the original neighbourhood U if necessary, we may write 

F,{u) = Fo{u) + eT{s,u) (2.5) 

with u G U C K and e € M'' by analogy with the original problem in Section [l] Note that FeiOx) 
is not necessarily zero when e ^ 0. 
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2.3. Uniform degeneracy 

It may happen that the map germ Fo\Kx is in some precise sense independent of a; e S": we shall 
make this more explicit in Section [3] below. In this case we say the normal degeneracy is uniform. 
For example, in the group orbit case with trivial isotropy (no nontrivial element of T keeps x € S 
fixed) the action of F provides a natural coordinate system for the whole of NS as a product bundle 
over S. In these coordinates the decompositions — (BL^ — Px®Rx are independent of a; £ 
and uniformity in a natural sense follows automatically. 

Machinery to study perturbations of functions and maps is provided by singularity theory. In 
the next section we introduce the key ideas that will be required for our purposes. 



3. Versal deformations 



For ease of exposition we make a simplifying assumption, namely that all of the bundles over S that 
are involved in the construction of the reduced vector field i^^ (including the tangent bundle TS are 
trivial bundles. This is automatically the case if S is contractible (a copy of M for example) and is 
usually the case in applications when 5' is a circle. Even when the assumption fails, we can apply 
the methods below to contractible open subsets of S and piece together the results subsequently. 
We remark later on how the methods would need to be adapted to apply intrinsically to nontrivial 
bundles. 



Let us therefore take K = S xM.^ and TS ® P ^ S xW^ x M.^ so that as in Section [2j] we 
regard f ^ as a family of smooth map germs 

parametrised by a; G and (small) e e M"^. Since we are here interested in the zero locus 
Zp = {(e, x, y) e M« X S* X M'^ : F,{x, y) = G 

for small |?;| rather than other geometric features of as a map, we are free to apply coordinate 
changes in a neighbourhood of Ox S'q = 0x5x0 that respect the geometry of the zero locus and 
the role of e as a parameter. A natural class of coordinate transformations to use is the class of 
contact equivalence,^^ or JC- equivalences. Singularity theory provides powerful tools for handling 
classification and geometry of map germs up to /C-equi valence, and the essence of this paper is the 
adaptation of those techniques to our class of bifurcation problem. It is in order to use the tools of 
singularity theory that we assume C°° smoothness of all data in the problem. Background material 
for the singularity theory in this section can be found in a range of texts such as [5j (6) i25j and from 
a slightly different viewpoint in [27] . 



3.1. Basic versal deformation 

We here introduce the main elements of the theory of /C-versal deformations of map germs that will 
be used in our analysis. To recall the concepts we shall use generic symbols x and F that are not 
the same as a; G S* and the map F^ of our main problem. 

1 Not to be confused with contact transformations from classical mechanics. 



Definition 3.1 Two smooth map germs f,g : M",0 W,0 are /C-cquivalcnt if there is a 

diffeomorphism germ (f) at the origin in M" and an invertihle p x p matrix A depending smoothly on 
xeR" such tha^ 

g{x)EEA{x)mx)). (3.1) 

Clearly the zero set for / is the image under (f) of the zero set for g . Since A is invertible any 
degeneracy of / at the origin is 'the same' as that of g . 

Definition 3.2 An r-parameter deformation of a map germ f : ]R",0 — > W,0 is map germ 
F -.W X M", (0,0) RP,0 such that F{0,x) = f{x). 

To keep notation simple we shall write _F as a map, the fact that it represents a map germ at (0, 0) 
being understood. 

Now for the key definition through which we classify types of bifurcation behaviour. 

Definition 3.3 Two deformations Fi,F2 : M.^ x M" ^ of f are /C -equivalent if there is a 
diffeomorphism germ $ at the origin in MJ" x M" of the form 

^a,x)^{i;{a),.p{a,x)) (3.2) 

and an invertible p x p matrix A varying smoothly with respect to (a, x) G M"" x M" such that 

F2{a,x) = A{a,x)Fio<^>{a,x). (3.3) 

Note that Fi, F2 are /C-equivalent as maps, but the coordinate change $ distinguishes between the 
roles of the variables x and parameters a . In particular, <I> takes the zero locus F^^iO) to Ff 1(0), 
and tp takes the parameter-space projection of F2^{0) to that of F^^{0). 

In the case when A{x) or A{a, x) is the identity matrix (hence can be deleted from the 
expressions (3.11 or (3.3 1 ) the maps or deformations are right equivalent {TZ- equivalent). It is easy 
to check that if two real-valued functions are 7?,-equivalent (up to a constant) then their gradient 
vector fields are /C-equivalent. 

Much of the power of singularity theory resides in the central notion of versal deformation: a 
typical map germ has an identifiable class of deformations from which all others can be obtained 
by suitable coordinate substitution. 

Definition 3.4 The r-parameter deformation F of f is /C-versal if given any q-parameter 
deformation G{c,x) of f , there is a map germ a : W,0 M'',0 with the property that G is 
IC-equivalent to the deformation 

a*F : (c, x) i-^ F{a{c),x). 

The deformation F is AI!-miniversal if there exist no IC-versal deformations of f with fewer 
parameters. 

In other words, the parameters Ui in the /C-versal deformation F can be expressed as suitable 
(smooth) functions ai of the given parameters Cj in order to recover G from F up to /C-equivalence. 
Note that in general a is not a diffeomorphism, since (apart from other considerations) q and r 
need not be equal. 

Remark. It is a consequence of the general theory that if F is /C-miniversal then F{a, •) and / are 
/C-equivalent germs if and only if a = S M'' . 

+ The symbol = indicates the two map germs are identical: it is not an equation for x. 
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If / is nonsingular (that is Df{0) has maximal rank) then the IFT imphes that every 
deformation of / is /C-versal. In what follows we therefore suppose that / is singular. 

Not every / necessarily has a /C-versal deformation, but the failure to do so is (in a precise 
sense) exceptional. Moreover, and most importantly, IC-veisal deformations when they exist can be 
chosen to take the explicit form 

F{a, x) = aiVi{x) + • • • + arVr{x) + f{x) (3-4) 

for a certain (nonunique) collection of polynomial map germs Vi : E" W , 1 < « < r, that can be 
constructed from / by a standard algebraic procedure. 



From Definitions 3.2 and 3.4 we deduce the following: 

Proposition 3.5 The zero set Zp := F~^{0) is a 'standard object' in MT" x M" associated with f : 
the zero locus Zq of any given deformation G of f as above is (up to difjeomorphism ^) the inverse 
image of Zp under the map (c, x) i~* {a{c),x). □ 

We are concerned with bifurcations, so the points of Zq of particular interest to us are those 
where the projection of Zq into parameter space fails to be a submersion: these are points where 
the local structure of the set of solutions x to G{c, x) — can change as c varies locally. 

Definition 3.6 The singular set of the deformation G : M'' x K" Rp is defined as 

Eg {{c,x) e M« X M" : rank D^G{c, x) < p}. 

The discriminant Aq is the projection of Zq n Y,q into the parameter space M'. 

Thus Aq is the set of c G M"^ for which there exists some x G M" with G'(c, x) = but DxG{c, x) 
has rank less than p. 

The importance of Y,q for us lies in the following key result much used in deformation 
(unfolding) theory. See e.g. (25, p. 49] or [17, Ch.ll,(18)] for a proof. 

Proposition 3.7 If G is a submersion, so Zq is a smooth submanifold o/M'' x M" of codimension 
p , then Y,Q is the set of singular points of the projection Zq <zW^ x M" M'^. □ 

Remark, li n < p (a case important to us) then T.q — Zq and the rank condition is vacuous. 



From Definition |3.4| it follows that if _F is a /C-miniversal deformation of / then, for any other 
deformation G of / , the discriminant is (up to diffeomorphism ip) the inverse image of the 
discriminant Ap under the parameter map a. The discriminant Aj? is also a 'standard object' and 
for polynomial / it is given by a finite set of polynomial equations in W . 

The policy statement below is the ideological basis for our geometrical treatment of bifurcation 
from S in this paper: 

The key to understanding the geometry of the zero locus of the deformation G of f is to 
understand the geometry of the associated map a from the given parameter space of G 
into that of a K-miniversal deformation F of /, and in particular the intersection of the 
image of a with the discriminant Ap . Generic assumptions about this intersection lead 
to conclusions about generic geometry of Aq. 

Because of the degeneracy along S we shall need to take extra care in implementing this policy. 
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3.2. Deformation along a manifold. 

There is a vital technical issue to deal with, before we can exploit the above ideas. In our context as 
described in Section [2T2| the full parameter space is M.'^ x S, with the e-deformation terms vanishing on 
0x5'. However, the notion of versal deformation as described above applies only in a neighbourhood 
of a single point in parameter space, whereas we would like to apply it in a consistent way over a 
whole neighbourhood of x 5 in M"? x 5. 

In standard proofs [6l [25] of the existence of /C-versal deformations the requisite coordinate 
changes are constructed by integrating suitably chosen vector fields, the existence of these vector 
fields being established by using the Preparation Theorem for smooth function germs. However, as 
pointed out in [HI (see Section 4.5, Example 4), there is a parametrised version of the Preparation 
Theorem which holds where the parameter belongs to a closed interval [0, 1] and where (as in our 
case of uniform normal degeneracy) the map germ to be unfolded is independent of the parameter. 
The proof of this theorem uses (as expected) a smooth partition of unity. The same proof applies 
equally well to a parameter belonging to any smooth compact manifold (there are no topological 
obstructions in view of the linear nature of the constructions). This enables us to work with the 
'global' parameter {e,x) g x 5 for |e| small. We can therefore state the following deformation 
result, customised for our purposes: 

Proposition 3.8 Let h ■.R'',0 ^ R'^+'', be a smooth map germ which has a JC-versal deformation 
H ■.M.'^xR'' ^ R'^+^. Let W be a neighbourhood of x Sq mR'^ x S xR^ and letF -.W ^ K'^+'= be 
a smooth map such that i^(0, x, y) = h{y) for all x G S. Then there exists a neighbourhood W C W 
of X Sq and a smooth map 

<^:W' ^R^ xR" : (e, x, y) ^ (V(£, a:), </>(£, x, y)), (3.5) 
as well as a nonsingular matrix A{e,x,y) of size d+k and defined on W' , such that 
(i) the identity 

F{e, X, y) EE A(e, x,y)Ho $(e, x, y) (3.6) 

holds for all (e, x, y) £ W ; 
(ii) the map R^ — *■ R^ : y t-^ y ^ (l){e,x,y) is a diffeomorphism for each (e, a;) E R"^ x S, with 

(/)(0, x,y) — y for all x G S ; 
(ill) V(0, x) = e M'' for allx e S . 

□ 



Putting together Proposition 3.8 and the expression (3.4 1 we now obtain the deformation result we 
require for application to the reduced bifurcation problem i^^ = as in (2.1 1. 

Corollary 3.9 Let the vector field Fq satisfy uniform normal degeneracy on S as in Section \2.S\ 
so that in coordinates {x,y) on the neig hbourhood U of So in K ^ S x M'^ we have Fo{x,y) = h{y) 



where F^ is the reduced map as in {2.1). Lf h has JC-versal deformation of the form 

H{a, y) = aivi{y) H arVr{y) + h{y) 

then, up to premultiplication by a nonsingular matrix depending smoothly on (e, x, y), we can express 
F^ in the form 

F,{x, y) = F{e, x, y) = ai(e, x)vi{y) + ■■■ + a^(e, x)vr{y) + h{y) e R'^+'' (3.7) 
where y t-^ y is a diffeomorphism germ at the origin in M*^ depending (smoothly) on {e,x) E R"^ x S 
and where also 

a,{0,x) = (3.8) 
for all X E S and i = 1, 2, . . . , r. □ 



11 



3.3. Global and local geometry of the discriminant. 

We continue to assume uniform degeneracy, and without loss of generality we take to be of the 
form (3.7) where after the coordinate change y y we drop the tilde on y. The identity (3.8 1 
implies that we can write 

a,;(e,a;) = h{x)s + Ci{x){s,e) + 0{\s\^) G M 

for i = 1,2, ...,r, where bi{x), Ci{x) are linear and bilinear maps — 
respectively. Thus 

b^{x)e ^^bij{x)ej and c,(x)(g, e) = cumix) £; 



(3.9) 



and 



Im 



where 



hj{x) 



(e,x)|e=o and cumix) 



dsidsm 



-(e,x)|e=o- 



(3.10) 
(3.11) 



We next write e in 'polar coordinates' as 

where IR+ = {t e M : t > 0} and S^^^ = {s e M"? : \s\ = 1}. From the smooth map 
a ■.R'^ X S : (e, x) i-^ (ai(e, x),. . . , ar{s, x)) 

satisfies 

ai{e, x) = pbi{x)s + p^Ci{x){s, s) + O(p^) G M 

for ? = 1, . . . , r. We study the geometry of the map 

b-.S"-^ xS^W : {s,x) ^ {bi{x)s,...,br{x)s). 

From this we obtain the geometry of the map a : M'' x 5 ^ M'' to first order in p as the cone on that 
of b . We anticipate that under suitable generic assumptions on b this cone structure will enable us 
to capture the true geometry of the discriminant 

of F near the origin. 

To investigate this closely we look at rays of the cone and their intersections with Ah as p 0. 
It will be convenient to do this using a 'blow-up' or (in reverse) pinching technique as follows. First 
define the pinching map 

T : M X ^ R"^ : (p, s) ^ (ps) 

and write 

T -.^T xid:Rx §9-1 X S -^R"^ X S, 

noting that since t{p, — s) = t(— p, s) we may restrict attention to p G M+ := {r G M : r > 0} if 
convenient. Next, let 



Dl :^ b'\A„) C X S. 

ont 
no( 

A^ ~ a-^{A„) cR+x §"-1 X S 



Generic assumptions on b will enable us to characterise the local structure of Dp which then 
determines local models for the structure of 
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where 

a:=aof :M+ x S^-^ x S^K'. 

(Note that here a plays the role of the map a in Definition 3.4 ) Finally, the application of f gives 
the local geometry of = a^^{Aff). 

The following elementary fact about t will be important. 

Proposition 3.10 Let 7 : M, ^ M x S'?"-'^, (0, xq) be (the germ of) a smooth path such that 

7(i) - ipt"' + 0(f"+i), so + tu+ 0{t^)) 

where p ^ G R and u ^ G TsgS"^^^ (so u J- sq in ). Then the image of S := t o ^ in 
is (close to the origin) the union of two arcs from the origin, each having contact of order 
(m + l)/m with their common tangent direction Sq . 

Proof Immediate, since T'^{t) = pt^so + O{f^+^)so + pf^+^u + 0(i™+2), so taking the 
parametrisation t' = t"^ gives 

U'\ J-' ,r{m+l)/m 

T^{t)—pSQt = pt U 

to lowest order in t'. □ 

Observe that if m is even then 5 has a cusp at the origin in W , while if m is odd then i5 is a 
1-manifold passing through the origin. 

The following geometrical construction will also play a key part in our investigations. 

Let H -.W y. M" ^ M^* be a /C-miniversal deformation of the singular germ h : M", — > M^, 0. 

Definition 3.11 A vector v gW is called a positive tangent vector to A/f at the origin if there 
exists to > and a differentiahle path 7 : [0,to) with 7(0) = G such that 

(i) j{t) G Ah for alltG [0,to) 
(n) 7'(0) = V. 

The positive tangent cone Th for Ah at the origin is the union of all the positive tangent vectors 
at the origin. 



We now turn to applying this machinery to the reduced bifurcation problem (2.1 1. 



3.4^. Branch points and bifurcation 

A& e gM3 moves away from the origin in a given direction we expect branches of solutions to — Q 
to emanate from certain points {x,Q) G Sq . 

Definition 3.12 A solution branch for F at x G S is a differentiate path z : [0, t^) — > Z p, C 

M« X 5* X M*^ for some to > such that 

(i) z'(0) ^ 

(li) z(0) = (0, X., 0) e X S'o 
(iii) z{t) e M« X S'o only ift^O. 

A solution branch z has direction s G S''^^ if the projection of z'{0) into is nonzero and in the 
direction s ; otherwise we say the branch z is tangent to S. If a solution branch exists at x G S 
with direction s then we say that x is a branch point with branch direction s . 
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Recall (see the Remark after Definition 3.6 1 that since k < d+ k the zero set Zp is here the same 
as the singular set T,p and projects to = a^^(A//) C M''. 

From the geometry of Ah we obtain a necessary condition for branching. 

Proposition 3.13 A necessary condition for xq Cz S to be a branch point with branch direction 
So G S''^^ is that 6(so,a;o) ™ the positive tangent cone Th for the discriminant Ah of the 
K,-versal deformation H of h . 

Proof. Let z be a solution branch at xq (z S with direction sq so that 
z{t) = {m,yit)) = ipit)s{t),xit),y{t)) eR'^xSxR'^ 



where p(0) = , p'{0) ^ and s(0) = sq , a;(0) = Xq. From (3.101 we have 



-a(z(t)\^^ = j^p(t)h{s(t\x(t)) + 0{p(tf 
and so 6(so,a;o) G Th since a(z{t)) e Ah for t e [0,to)- 



p'(0)6(so,a;o) 



(3.12) 

□ 



In the examples considered below the set Th turns out to be either a hyperplane or a half- 
hyperplane. Generically we expect the map b : S'~^ x S* ^ R'' to be transverse to Th in an 
appropriate sense, and the transversality allows us to establish the existence of branch points and 
to describe the local solution structure regardless of the terms 0(e^) in a(e,x). 

Because of the technical complications of the relevant singularity theory, the approach described 
above is useful in practice only for small values of d and/or k. However, the case A: = 1 is most 
significant, as it represents the least departure from normal nondegeneracy of S* , while the case 
d = \ naturally occurs in the study of periodic orbits of differential equations where there is an 
automatic circle group action on the space of periodic functions (see Section |6]) . Therefore the 
simplest but most important case d = k = \ merits detailed investigation. 



4. The fundamental case d = k=\. 



Here we have a 1-dimensional equilibrium manifold S of constant corank 1. Throughout this section 
we continue to assume that the normal degeneracy is uniform; some branching analysis applying to 
the general case is given below in Section [Sj 



To simplify notation we drop the tilde throughout and write (3.7 1 as 

F,{x, y) = F{s, X, y) = g{e, x, y) + h{y) e (4 1) 

with multiparameter e €W and m = (x, y) € S'xM since here k = 1 and d+k = 2. Each component 
of h{y) is without constant or linear terms because of the normal degeneracy of the zero set Sq. 
Supposing that at least one component has a nonzero Taylor series at 0, it is a straightforward 
exercise (see [2S] for example) to check that h is /C-equivalent to the germ y i-^ (0, y™) for some 
integer m > 0. We first take the least degenerate case m = 2 and study this example in detail 
in the geometric framework outlined above. We then show how this approach extends to general 
m and illustrate in particular how it recovers and sets in a geometric context the analysis of this 
problem for g = 2 given by Hale and Taboas [25] . 
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4-1- The case m — 2 : quadratic degeneracy. 
Here 

F{e, X, y) = g{e, x, y) + h{y) = .g(e, x, y) + (0, y^) G 
for (x, y) e 5 X M with dim 5 = 1 , and e e M?. 

4-.1.1. Geometry of the zero set. A /C-miniversal deformation of h{y) = (0,?;^) is 

H{a, y) = H{ai,a2,a3,y) = (ai + a2y, as + y^) 

(see [B],[1S]) and so by Corollary |3.9| up to /C-equi valence we can suppose F to have the form 

F{.£,x,y) = {ai{e,x) + a2{e,x)y ,az{£,x) + y'^) (4.2) 



for \y\, \e\ sufficiently small and for all x ^ S, where a(0, x) = G R for all a; G S*. As in Section 3.3 



in order to describe the bifurcation behaviour of F we need to understand the geometry of the 
discriminant Aj? = a^^(A/f) C M"^ x 5 where (as is easily seen by eliminating y from H{a,y) — 0) 
the discriminant A/f is given explicitly by 

af + alas = 0. (4.3) 

This is the equation of a Whitney umbrella in M^: see Figure [2] 




Figure 2. Whitney umbrella: locus a| + a| 03 = showing the stratification {To, Ti , Tj , T2}. 



From Section |3.3| the key intermediate step is to study the geometry of 
Dl := b-\AH) e S"^"^ X S 
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is as in (3.11 ). Following Proposition 3.13 we first identify the set b ^(T) 



where b : §«-i x S - 

where T = Th is the positive tangent cone to Ah at the origin. In the present case 



T={ae 



ii = 0,a3<0} 



since any smooth path in A^f from the origin can be checked to have contact of order at least | 
with the plane ai — 0. 

We now make the following generic hypothesis, which is at the heart of our whole approach to 
bifurcating from a normally degenerate manifold: 



(HI) The map b : S^-^ x S 



is transverse to T. 



The notion of transversality to a variety such as T rather than to a smooth manifold assumes 
a stratification of the variety, that is a suitable decomposition of the variety into finitely many 
manifolds of varying dimension, the smaller comprising the boundary of the larger. See Ch.3,§l] 
for example. In this case there is an obvious and natural stratification {Tq,Ti^T2} where Tq is the 
origin, Ti is the a2-axis with 02 9^ , and T2 is the part of the (02, a3)-plane with 03 < : here 
dimTi = i for i — 0,1,2. However, we see from Figure [2] that since our ultimate interest is in 
geometry of Aj? we should also take account of the fact that the negative os-axis T{ is a natural 
1-dimensional stratum in Ajj. We hence include this half-line T[ as a 1-dimensional stratum for T 
when interpreting the hypothesis (HI). 

Since Ah and T are 2-dimensional objects in M.^ the essential geometry of transversality is 
captured by taking q = 1 or 2 so that dim(§^~^ x S") = 1 or 2. We first study carefully the case 
q — 2, and comment on the cases q = 1 and q > 3 later. 



4.. 1.2. The 'typical' case q = 2. We take S compact for ease of exposition; adjustments to the 
noncompact case are straightforward. 

The transversality hypothesis (HI) implies that b^^{T) consists of a smooth 1-manifold with 
boundary, thus a finite collection of closed arcs or embedded circles in x S. The set of end 
points of these arcs is b~^(Ti), while interior points constitute b^^(T2) with possible intersection 
points 6~^(r{). Note from the geometry of Ah and linearity of b{s,x) in s that if {s,x) is an end 
point then so is (— s, x). Moreover, the intersection of 5*^ x {x} with Ap lies in the semicircular arc 
{s G : &3(s, x) < 0} since all of Ah lies in the half-space 03 < 0. A schematic representation of 
this geometry is shown in Figure |3] 

For {s,x) in a local coordinate chart around {sq,xo) G x write {s,x) = (sq + v,xo + w) 
where {v,w) e M x M. Using (3.101 we have 

a(p, v,w) ^ ao f (p, s, x) = pb{v, w) + p^c{v, w) + 0{p^) e (4.4) 

where now b{v,w) is shorthand for b{s,x) and analogously c(w, w) stands for c{x){s,s). From ( |4.3[ ) 
we see a{p, v, w) G Ah when 

{pbi + p'c^f + {pb2 + p'c2)\ph + P^C3) + 0{p^) = (4.5) 

which after cancelling p^ reduces to 

(3{p, V, w) bj + (26ici + blb3)p + d{p, v, w)p'^ ^ (4.6) 

where 6^, q are components of b(v^ w), c{v, w) and where d{p, v, w) is a smooth real- valued function. 
Now we study three possibilities, illustrated in Figure |4] 
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S 



Figure 3. Schematic representation of b~^{T) (dashed curve) and Dp = b~^{A[j) (solid curve) 
in the case d = k = l, n = q = 2 (Whitney umbrella) showing end points (black dots) and 
intersection points (white dots) . The whole of Dp lies in a strip of width 7r between two parallel 
curves and has quadratic tangency with those at the end points. The 1-dimensional manifolds S 
and §^ are each represented as closed intervals. 

(i) {so,xq) is an interior point: &i = and 62^3 7^ 0. 

The IFT at {v,w) = (0,0) implies that ii p := h\ ~ (3 then the map (pjVjW) 1-^ {p,v,w) is 
a local difTeomorphism at the origin. In {p,v,w) coordinates the locus f3 — is given hy p — bi, 
which (since (0,0) is a regular point of the function 61 by hypothesis (HI) of transversality to T2) 
is a smooth surface in M.^ quadratically tangent to the curve &i(w, w) = 0. 

Finally, to reconstruct the geometry in the original (e, x) deformation coordinates, we apply 
the pinching map f which takes the surface /3 = locally to a surface with a cusped ridge along 
0x5; the tangent direction of the cusp at (0, xq) is sq G S^, while as x moves away from xq in S 
the tangent direction s of the cusp is given by the smooth curve 

Bi = {{x, s) eS^ X S : biix, s) = 0} 

through {so,xo) in §^ x S. See Figure [i] 

(ii) {sq,xo) is an end point (denoted (sq,Xq)): 61 = 63 = 0,62 7^ 0. 

From the hypothesis (HI) of transversality to Ti we have '^g'^^ ''^^ 7^ at (0, 0), and so the map 
^ : {v,w) {i,fv) (5i,2ci&i + bjbs) 
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Figure 4. Schematic diagram of the discriminant Ap with cusped ridge along 0x5, showing 
the structure close to xo,Xq,Xq on S determined by an interior point (so,xo), end point {s'i^,x'q) 
and intersection point {sq,x'q) (indicated simply by so etc.) respectively. 




s 

Figure 5. The projection fi : Ap — > close to (0,xi) S M'^ X S where (si,a;i) is an interior fold 
point of the curve Bi . The cusp curve Ci projects to a cusp curve Ci (the bifurcation set) in 
tangent to the si direction. The numbers 0, 2, 4 indicate the number of solutions {x, y) a K to 
F{e,x,y) =0 close to (xi, 0) £ K for e close to the origin in R^. 
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is a diffeomorphism germ at the origin. In new coordinates the function (3 takes the form 
/3(p, V, w) — Q{p, V, w) + 0{p^) where 



+ wp + ep^ 



(4.7) 



which is a nondegenerate quadratic form regardless of the value of the coefficient e . The locus 
Q = is a cone C tangent to the (-D, i(;)-plane along the li-axis. By the Morse Lemma there is a 
local diffeo morp hism (with derivative the identity) that takes the locus /3~^(0) to the cone C and, 
moreover, (4.7 1 shows that j3~'^{Q) is quadratically tangent to the (w, i())-plane along the smooth 
curve V — hi ~ Q. Finally, by Proposition 3.10 the pinching map f takes /3~^(0) to a surface with 



|-power cusped ridge along x S* and terminating at (0,Xo), reminiscent of a ship's prow: see 
Figure |4j 




Figure 6. The projection fi : Ap 



close to (0, Xg) £M.'^ X S wliere {sq,Xq) is an end point 



of the curve Bi . The curve Cj projects to a pair of arcs Cj (the bifurcation set) in 'M? tangent to 
the s'q direction. The numbers 0,2 indicate the number of solutions {x,y) K to F{e,x,y) = 
close to (xq, 0) £ -fC for e close to the origin in M.'^. 



(Hi) (sqj xq) is an intersection point (denoted (sq, a;Q)): 6i = 62 = , 63 < . 

By the hypothesis (HI) of transversality to T[ the map {v,w) (^1,62) is a local 



diffeomorphism at (0,0). From 
(fli + ia2){ai - 

where £ = 



ia2) 



3| with 03 < we can express Ah as 


with £ = 



■^/plfoaj. Hence for fixed 



/1 03 1, which on cancelling p~ becomes bi ± £b2 
p > the expression (4.6 1 represents a pair of smooth curves crossing transversely at the origin, 
with tangent directions both tending to that of bi{v, w) = as p —> 0. See Figure [i] 
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4-. 1.3. Bifurcation geometry when q = 2. Wc have described above the three generic types of local 
geometrical structure found in the discriminant Ap , where we recall that (since I = d < p = 
d + k — 2) the discriminant is the projection into x S* of the zero set: 

Ap = {{e,x) £R^ X S -.By eR, F{e,x,y) = 0}. 



However, it is important to keep sight of the fact that in terms of the original problem (2.1 1,(4.1 1 it 
is (a;, y) that are the variables while e is the perturbing parameter. Therefore to find the bifurcation 
behaviour of solutions (x^y) as e varies in we must also study the geometry of the projection 
p : Ap —> into the e-plane and in particular its singularities. The image of the singular points 
(i.e. the set of singular values) is the bifurcation set; this can also be interpreted as the apparent 
outline (or apparent contour: see [H].|22j) of Ap when viewed along the 1-manifold S. 

The analysis of the previous section showed that Ap has a cusped ridge along x S* with 
the direction of the axis of the cusp (indicated on the unit circle C M^) given by the curve 
Bi = b^^{0) C §^ X 5*. (For any interior point x £ S there can be only one s e §^ with (s, x) £ Bi, 
so cusps may impinge on (0,x) £ x S from only one direction, while at end points they do so 
from two opposite directions.) Therefore the singular points of /i|_Bi : i3i — > S"'^ play a key role in 
the structure of the bifurcation set for F. Accordingly, we make the following generic hypothesis: 

(H2) The map /i|-Bi : Bi ^ E>^ has only fold singularities, and these are at interior points of Bi . 

A fold singularity of is characterised in local coordinates {v, w) on §^ x S* by 

1^ = 0' -Q^^^ ^* («,«;) = (0,0). (4.8) 

Following the classification in the previous section we now consider in turn the three types of point 
on the curve Bi. 

(i) Interior points. 

The first result confirms what is evident from a simple sketch. See Figure [5j 

Proposition 4.1 Let {si,xi) he an interior point of Bi C x S" such that the projection 
p\Bi : Si ^ §^ has a fold singularity at (si, Xi). Then the apparent outline in R^ of a neighbourhood 
of (0, Xi) in Ap viewed along S is a curve Ci with a cusp of order | at the origin and with tangent 
direction si there. The number of solutions {x,y) to Fg(x,y) — increases by two as e crosses each 
arc of Ci in the appropriate (and same) direction of rotation about the origin. 

Proof. In local coordinates (p, v, w) the apparent outline Fi of A\, viewed in the direction of S is 
by definition given by 

Fi = |(p,w) £R+xR:3w £ S , P{p,v,w) = ^{p,v,w) = 

ow 



and we have already seen in Section 4.1.21 (?) that in suitable coordinates (in which the role of w 

" l^=26if^ 

aw ^ aw 

dw 



is unaffected) we may take [3 — b\ — p. Hence 5^ = 2bi ^ which vanishes where bi — (so then 



p = if /? = 0) and also where |^ = 0. By the IFT and the fold condition (|4.8|, the latter locus 



is R+ X B'l where B'l is a smooth curve through the origin in not tangent to the w-axis. The 
intersection Fi ofM+ x B'l with Ap is a smooth curve quadratically tangent to the (u, w) plane. 



By Proposition 3.10 the curve Ci := f (Fi) has a ^-power cusp at (0, xq) and projects under ^ to a 
curve Ci with the properties described. □ 
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Note: The geometry of the curve Ci in = a^^(A^f) and the expression (4.2 1 show that as 
£ crosses the arcs of Ci one pair of solutions (a;, y) is created (or annihilated) close to y — yi where 
yi = ■\/— 03(6, x), and the other pair close to —yi- 

(ii) End points. 
The analysis in Section 



4.1.2 



I showed that at an end-point the surface Ap = (3~^{0) has a Morse 
singularity, thus a local cone structure, with the cone tangent to the (u, z/;)-plane along the curve 
Bi. Let X = f§ • From (4.7 1 we see that x — when p = and 61 = , that is on the curve Bi. 



Differentiating the expression (4.6 1 with respect to w we then observe that provided the generic 
condition 



(H2') ^(26ici 



qbs) + 



holds at (0, 0) then §^(0, 0) 7^ and so by the IFT the locus x ^(0) locally has the form of a graph 



Now 

aw 



9x / 9x 



I -J^ which is nonzero at an end point because on 



dw I dp 



of a smooth function p 

Bi we have || = 2(|^)^ 7^ by the hypothesis (H2). Thus the smooth surface x^^O) C M x 
cuts the (w,u')-plane and hence the 'cone' Ap transversely along Bi (away from the origin) and so 
must cut the cone again along another smooth curve T'l through the origin and transv erse t o the 

a Ci 



3.10 



w,w-plane. Then C[ = f(T'i) is a smooth curve whose image under fi is by Proposition 
arc C'l through the origin in M^. See Figure [g] To summarise: 

Proposition 4.2 Let (sq,Xq) be an end point of Bi CZ x S and assume that (H2' ) holds. Then 
the apparent outline in of a neighbourhood of {Q, x'^) in Ap viewed along S is a arc C[ through 
the origin, across each branch of which the number of solutions {x,y) to F^{x,y) — increases by 
two, taken in the appropriate (and opposite) directions of rotation about the origin. □ 

Remark. The arc C( is typically not the same as the 'end-point' curve a~^(Ti) which itself 
corresponds only to solutions with j/ = 0. Thus soon after a pair of solutions is created across 
C'l one of them will pass through y — Q, that is through 5*0 . However, we are here solving the 



reduced equation (2.1 1, so the corresponding equilibria of ( 1.1 1 for e ^ need not lie on S although 



and the zero set of 61 is the 



they will have zero if-component and therefore lie in L. 
(Hi) Intersection points. 

Intersection points are zeros of the map (61,62) : x 5* ^ 
curve Bi . We make a further generic hypothesis: 

(H3) The function 62 does not vanish at the singular points of fJ.\Bi . 

An intersection point (sq, Xq) corresponds to an arc of self-intersection of Ap emanating from 
{0,xo) : see Figure [ij This arc projects under p, to an arc C" from the origin in K^, and (so,a;Q) 
is a regular point of p\Bi by the hypothesis (H3). Thus no bifurcations occur across Cf , but two 
solutions of the reduced equation of opposite sign in y pass through y ^ (that is, through Sq ) 
simultaneously as e crosses C", so the corresponding solutions of (1.1 1 pass through L. 



We have now studied the three generic types of local bifurcation behaviour, and we assemble 
them into a global statement. Recall that d — k = 1 (so the compact manifold S is 1-dimensional 
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and is degenerate with corank 1) and the degree m of degeneracy is 2 . The perturbation parameter 
£ Hes in M^. 

Let K denote the local involution in a neighbourhood of S in M that in terms of normal bundle 
coordinates TS ® K ® L corresponds to (x, y, z) ^ {x, — y, z). 



Theorem 4.3 Under the generic hypotheses (H1),(H2),(H2'),(H3) the bifurcation set for (1.1) 
consists of a finite set of arcs emanating from the origin in . The arcs are of three types: 
fold, end and intersection arcs. These correspond to the following bifurcation behaviour as e rotates 
about the origin in M? in (locally) the appropriate direction. 

Interior arcs. These occur in pairs with ^-power cups at the origin. As e crosses the first arc 
a pair of solutions is created away from S by a saddle-node (fold) bifurcation, and as e crosses the 
second arc a second saddle-node pair of solutions is created close to the K-images of the first pair. 

End arcs. These occur in opposite pairs forming curves through the origin, and each arc 
is accompanied by an arc of a~^{Ti) with the same tangent at the origin. As e crosses an end arc 
a pair of solutions is created away from S in a saddle-node bifurcation; as e crosses the companion 
arc one of these solutions passes through L. This behaviour is repeated at the opposing end arc and 
with the same sense of rotation about the origin. 

Intersection arcs. As e crosses an intersection arc with nonzero speed two solutions of the 
reduced equation with K-components of opposite signs pass through L at the same point and with 
nonzero 



This theorem in particular recovers and provides a geometric context for the results of Hale and 
Taboas 29J for the case m = 2. It also extends those results by taking account of self-intersection 
points of the determining curve Bi. See Figure [t] and compare Figures 1 and 2 in |29j . 

4.1.4- Cases q ^ 2. For g = 1 (that is e e M) the transversality hypothesis (HI) implies that the 
image of 6 : x 5 ^ M'^ meets T only at points of T2 and that is a regular value of 61 : S" ^ M. 
Hence 6^^(0) is a discrete set of points {(±1,0;^)}. (This may be visualised in Figure|3]by taking a 
pair of vertical lines corresponding to fixed values ±s € S^.) To each Xi with either e > or e < 
there corresponds a pair of solution branches in Ap cRx S emanating from (0,0;^); the branches 
are such that locally x — Xi varies as |e|2. In particular this describes the branching behaviour 
for q = 2 when e leaves the origin along an arc not in a direction tangent to an interior, end or 
intersection arc as in Theorem 14.31 

For (7 = 3 some interesting new geometry arises. For x E S the map a{-,x) : M.'^ ^ M.^ is 
typically nonsingular at the origin, although generically there is a set of isolated points {xj} C S 
where the derivative Dj of a{-,Xj) has rank 2. For x ^ {xj} and for sufficiently small |e| the 
discriminant Ap = a^^An intersects M.^ x {x} in diffeomorphic copy {A|f} of Ah (more precisely, 
a neighbourhood of the origin in Ah , which amounts to the same thing). To capture the structure 



of the bifurcation set we study the map 6 : §^ x S" ^ M as in (3.9 1 and the geometry of b ^(T) 
where T = Th . Writing b^{s) = b{s,x) we note that b^T) = B^yj Bl where B° is a 

semicircular (open) arc b~^{T2) of a great circle Cx in with antipodal end-points B^ = b~^{Ti), 
while Dp{x) := 6~^(A//) is a figure eight curve on passing through these end-points and with 
self-intersection at bx^{T[). As p ^ the intersection pS^ n A^ flattens to pb^^{T) at rate . 

In this setting the bifurcation set in M.^, which is the set of critical values of the projection 
p : Ap C S'^ X S —> M.^ , is conveniently interpreted as (to first order in p) a cone on the envelope 
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Figure 7. Schematic bifurcation diagram for d = k = l, q = 2 and m = 2 (see text), showing 
end arcs (unbroken lines) and fold arcs (long dashes) through the origin in the parameter space 
R^. Intersection arcs (short dashes) are also shown. Numbers indicate the number of solutions in 
each region. Note that for each opposite pair of fold arcs the two folds approach from the same 
side in the plane. 

E of the family of curves {Dp{x)^ ^^^'m The curve E separates points s S according to the 
number of points x E S with s G Dp (x) . The generic local geometry of envelopes of smooth 
curves in (or S^) is well known (see ^ for example): the envelope consists of smooth fold curves 
across which Us changes by 2, meeting at isolated cusp points. These local configurations are stable, 
so that the same description applies to the envelope of the family {pE>^ n A^^^^^ (where higher 
order terms in p are not neglected) for p nonzero and sufficiently small. 

We could leave the discussion at this point, but further useful insight is obtained by taking 
account of the role of b~^(T), which we regard (compare Section |4| as a 'first approximation' to 
Dp{x). Let C be the envelope of the great circles {Cx}x£S and C C be the envelope of the 
semi-circular (open) arcs Let be the locus of end-points {B^} for x £ S. It is 

straightforward to check that if sq G C fl and is a fold point of C then B^ is a smooth curve 
quadratically tangent to C at Sq where an arc of terminates: see Figure ^a). When we now 
replace B° by the figure-eight Dp{x) then the 'approximate' envelope i?° U B^ becomes a generic 
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Figure 8. Envelopes in sketclied in the plane: (a) Local configuration of tfie envelope ij" 
of arcs _BJ! and the locus of their endpoints, (b) corresponding configuration of envelope of 
smooth curves Dp{x). 



envelope configuration E consisting of a fold curve close to B'^ which continues as a fold curve 
close to one branch oi B^, together with another fold curve close to 5° which meets a fold curve 
E^ close to the other branch of B^ at a cusp point close to sq : see Figure [sj^b). The numbers 0, 2, 4 
indicate relative numbers for s in each of the three complementary regions of E locally. Finally, 
local stability implies that this configuration occurs in pS^ fl Ai;- for p sufficiently small. 

Features of the earlier q = 2 case can be seen by intersecting the above with a typical great 
circle we find interior points in pairs (corresponding to the pair of branches of E close to Co) 
and isolated end points (corresponding to the branches of E close to Ci ). The q = 3 configuration 
shows how in the presence of a further bifurcation parameter these features generically coalesce. 

Returning to the exceptional points {xj} on S where lankDj = 2, let us make the generic 
assumption that for each Xj neither 1-dimensional stratum Ti or T{ of T lies in the range of Dj. 
It follows that for a; € 5* close to xj the image fe^rfS^) is an oblate ellipsoid with 'equatorial' plane 
close to the range of Dj. From this and Figure^ it can be seen that Dp{x) has the figure-eight 
configuration with one loop much larger than the other. As x passes through Xj the smaller loop 
vanishes (the end point coalesces with the intersection point) and is re-created at the antipodal end 
point; the transition curve D],{xj) in has a |-power cusp at both its antipodal end points. 

Provided Xj does not correspond to a cusp point of the envelope E the local fold structure of the 
envelope itself (that is, the bifurcation set in §^) is unaffected as x passes through xj , although the 
location of corresponding solutions to the bifurcation problem is controlled by this extra twist. The 
bifurcation structure at Xj corresponds to one of the generic singularities of 1-parameter families of 
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plane curves at the envelope: see 

For <7 > 4 generically for all a; G 5 the map a{-,x) : M'' ^ M"^ is a submersion at the origin: thus 
the local structure of Aj- is a product M.'^~^ x Ap where Ap is as in the case q = 3. Nevertheless, the 
envelope of the family {Dp{x)} in S*~^ may have additional singular structure, corresponding 
to higher-dimensional generalisations of Figure |8] which we do not investigate here. 

We now turn to cases of greater degeneracy degree m > 2. The essential geometry, while in 
principle more complicated, has much in common with the case m — 2 above. 



4-.2. The case m > 2. 
Here 

F(0,x,2;) = %) = (0,y'") 

and this germ has /C-codimension 2m — 1 (see [5] , |25j ) with AI!-miniversal deformation (using more 
systematic notation than previously) of the form 

H{a,y) = {ao + aiy -\ ha 

where a = (a, a) = (oq, . . . , a^.i, ao, . . . , 0^-2) G 
where Rm is the resultant of the two components of H{a, y) and which has the following expression 
as a determinant that for clarity we illustrate for m — 4, the general pattern being analogous: 



-12/" \ ao + aiV + 



.,"1-2 



p2m— 1 



am-2y - - + yn (4.9) 

The discriminant A/^^ is given by i?„j = 
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Here the blanks represent zero entries. 

Proposition 4.4 The resultant Rm has the form 

R,n = a" + (^oPrn + (-l)"ao(a5" + aiQ,n) + 00(4" + Rm) 

where Pm,Qm o,re polynomials each of least degree m and Rm is a polynomial of 
m + 1. 



(4.10) 



(4.11) 

at least 



Here the polynomials are in the variables ai , aj and least degree means the degree of the lowest-order 
term. 

Proof. The lowest-order term is a™ from the main diagonal. Putting ao = we then have 
Rm = {—l)™'aoRm-i where R'm_i is a determinant of order 2m — 2 illustrated in the case m = 4 
as follows: 

ai a2 as 

ai a2 aa 

ai a2 as 

ai a2 as 
ao ai a2 1 

ao ai a2 1 



(4.12) 
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The lowest-order term in R^-i is a™, and then if ai = the ao hi the first column is multiplied 
by a determinant whose lowest-order term is (— l)'"a™. □ 

Since the term of lowest order in i?,„ is a™ we see immediately: 

Corollary 4.5 The positive tangent cone Th lies in the hyperplane ao = 0. □ 

As in Section |3.3| we take polar coordinates e — ps and write 

(a, a) — {a{ps, x),a{ps, x)) = p{b{s, x),b{s, x)) + p'^{c{s, x), c(s, x)) + 0{p^) 



where {b,b) : S^-^ x 5* ^ M^m-i ^^^^ likewise (c,c). We then find from ( |4.1l[ ) that i?„ = p"^Rm 
where 

R,n = 6, 6) = 6^ + p((-l)™&o&r + ^o^^m) + + 0(p3) (4.13) 

in which P,„, Qm+2 are homogeneous polynomials of degree m, m -I- 2 respectively in the variables 
{bi, bj). Following the same procedure as in the case m = 2 we first study the geometry of the set 

A), {{p,s,x) e M X X S : Rm{p,b{s,x),b{s,x)) = 0} 

since Rm{p, b,b) — Q ior p ^ Q precisely when (a, a) G A// , and then apply the pinching map f to 
obtain = ^(A),) x S. 

We shall require some transversality hypotheses analogous to (HI). Rather than formulate a 
general statement (which we do not need and in any case cannot carry out without specifying a 
stratification of Th) we state the appropriate hypotheses as we come to them; they mirror those of 
Section I4TT1 

Suppose (so, xq) satisfies 6o(so: xq) = 0, and take a local coordinate chart (s, x) — (so4-w, xq+w) 
on §'~^ X 5 as before. As in the m — 2 analysis we look at three cases. 

(i) bobi 7^ at {sn,xo). 

Here by the IFT the locus R„i = locally has the form 

p = kob"^ + Un+i , fco ^ 0, (4.14) 

where 6o = &o('So + v^xq + w) and where Im+i denotes terms of degree at least m -I- 1 in w, w. We 
now make the generic hypothesis 

(Hl„i(i)) (so,a;o) is a regular point of the function bo- 

Here Ap is locally a smooth hypersurface in M x S"^^^ x S having contact of order m with 
X S'^^ X S along the codimension-1 submanifold 

Bn {{s,x) G S"-^ x 5 : bo{s,x) = 0}. 



Applying the pinching map r, remembering that t{p,—s) = T{—p,s) and using Proposition 3.10 



4.1.2 



we obtain Ap locally as a C"'^ hypersurface in K'^ x 5' containing x and having contact of order 
(m+l) /m along Ox S with the 'direction locus' in M'' x 5 determined by Bq . Thus for g = 2 , when 
m — 2 the discriminant A^^ has a |-power cusped ridge along x S" as already seen in Section 
while for m = 3 it is a submanifold of x passing through x S*. 

(a) 6o = 0,6i 7^ at {so,xo). 
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Here we make the next generic hypothesis: 

(H1„,(m)) {so,xo) is a regular zero of the map (5o,&o) ■ ^''^^ x 5 — + M^. 

First take q = 2. The transversahty hypothesis (Hlm,{ii)) means that we can regard the pair 
of functions (5o,6o) as local coordinates on x 5 close to (so,a;o). We can then write Rm in the 
form 

p"' + p{p + cp + 0{2)) (4.15) 

where {p,p) '■— {bo, (— &i)™6o + boPm) and c is a constant. The quadratic terms are nondegenerate 
regardless of c, so by the Splitting Lemma (see [I7j,|22 §4. 5], [33]) we can suppose the coordinates 
are such that Ap has the form 

p"' + p{p + cp) = 0. (4.16) 

The description of this locus depends on whether m is even or odd. 

The case when m is even is relatively easy, li m — 21 then we write t — p^ and the left 



hand side of (4.16 1 becomes a nondegenerate quadratic form in {p,t,p): the zero set is a cone with 
vertex at the origin, tangent to the (t,p)-plane and with the p-axis as a generator. On reverting to 
(p,p,p)-coordinates the cone retains these features, although the contact with the plane p = is 
now of order m rather than quadratic. Applying f then gives Ap with a 'ship's prow' geometry at 
(0, xo) in R'^ X 5 as in Section 4.1.2| (see Figure |4| but with cusp contact of order (m + l)/r 



im. 



The case when m is odd is a little harder to describe pictorially. Some insight is gained by 



regarding (4.16) as the equation of a curve in the (p,p)-plane for each fixed p. As p increases 
through zero, the order-m tangency of the locus with the plane p = along the p-axis reverses its 
orientation as a 'bulge' passes through the origin from positive to negative p where as before we 
interpret (— p)s as p(— s) when p < 0. A 3-dimensional representation is attempted in Figure [9j 

For the general case when q> 2 the local structure of A^? is the cartesian product of the above 

with m.i-'^. 

{Hi) 6o ^ 0, 6i = at (sq, 2^o)- 



Substituting ao — ai =Q into the expression (4.11 ) for i?,„ we obtain Rm = a^Dm where Dm 
is a polynomial whose term of lowest order is a™. Therefore with q = 2 we make the further generic 
hypothesis 

(Hl™(Mi)) 62(50,2:0)7^0. 

It then follows that (e, x) ^ for sufficiently small ja; — xq] and \e\. Thus there is generically 
no analogue of intersection points in Aj? when m > 2. 

When q > 2 there exist generically points at which bo — bi — b2 — 0, which we do not 
investigate further here. 
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Figure 9. The solution locus Aj^ given by Rm = for m = 3, g = 2 at a point where 
i>0 = bo = 0. Here the picture is constructed by taking i?3 = + p p + where p,p are the 
coordinates p = 6o> P = (— l)"*?)™ ^0 + boPm (see text). 



4-.2.1. The bifurcation set for m > 2. Having found the local form of the discriminant A^? cM'xS' 
under generic assumptions, we are now in a position to describe the generic structure of the 
bifurcation set, that is the ap parent outline of Ap viewed along S. Recall that Ap is the locus 
Rm = where Rm is given by (4.131 with {a, a) = p{b,b). Again we work first with q = 2. 



Taking local coordinates (s, x) — (sq + v,xo + w) on §^ x we seek common solutions to the 
equations 

d 



Rm — , 



dw 



Rm — 0. 



As before, we consider two cases. 



(i) 60^1 7^ at (50,2^0)- 



Here the argument parallels that of Section 4.1.3 We make the generic hypothesis 

(H2m(«)) The projection /i|i3o : -Bo ^ has only fold singularities. 

If ( So : 2:0) is a regular point of /i|i?o then no local bifurcations occur for any solution branches 
in directions sufficiently close to sq • If (sqj xq) is a fold point of /J,|i?o then from (4.14) and following 
Section 4.1.3 we find that the apparent outline of in the direction of S is locally a smooth curve 
To through (0, so) and having contact of order m with the (u, ti;)-plane, so that in particular when 
m is odd it has only one branch for p > Q while for m even it has two. Applying the pinching map 
f to A^p and projection fi gives a curve Cq through the origin in parameter space having contact 
of order (m + 1) /m with its tangent direction sq • The number of solutions (x, y) changes by two 
as each branch of Cq is crossed. 
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(m) 6o = 0, 6i 7^ at {sq^xq). 

Here we make the further generic hypothesis 

(H2™(zz)) ^((-6i)™^ + P„|^) ^Oat (so,xo). 

In terms of coordinates {p,p) used above, this simply states that 6 ■= ff and ?7 := §f are both 
nonzero at {p,p) = (0,0). From (4.15) and again disregarding the higher order terms as in (4.161 
we have _ 

dR 



dx 



so that ^ 



: Ap™ ^ with constant A 
we obtain 



-mC/?7 7^ at (0, 0). Substituting 



is solved as p = p(p) 
into ( |4.16| and deleting the factor p" 

p+ Ap + 0(p"~i) = (4.17) 

which represents a smooth curve B'q through the origin in the (p,p)-plane, transverse to both axes 
and not tangent to (^, rj) since m^l. Projecting the graph of p\B'q to a plane orthogonal to (^, ry, 0) 
then gives the apparent outline of viewed along S : this curve Cq has one branch for p > 
and one for p < when m is even and two branches for p > when m is odd. After applying the 
pinching map r the curve Cq becomes a bifurcation curve Cq having contact of order (m + 1)/to 
with its tangent line at the origin. For m even the number of solutions changes locally from to 2 
as the curve is crossed; for m odd it changes from 1 to 3 and back again in a hysteresis loop as the 
two branches are successively crossed. See Figure [TO] 




Figure 10. The discriminant for m = 3 viewed in a typical direction not parallel to p,p 
axes: apparent outline Cq (bifurcation curve) in becomes a cusped curve at the origin after 
applying the pinching map t. 



We summarise the bifurcation results in the following theorem. Recall that we are in the same 
dimensional setting as Theorem |4.3| but now with degeneracy degree m > 2. 
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Theorem 4.6 Under the generic hypotheses (lll„i{i) — (Hi)) anrf (H2m(i) — (m)) the bifurcation set 



for (1.1) consists of a finite set of arcs emanating from the origin inM?, each making contact 
of order ^^^Sl^I yjUli Hg tangent ray at the origin. The arcs are of two types: fold and end arcs, 
corresponding to the following bifurcation behaviour as e rotates about the origin in in (locally) 
the appropriate direction: 

Fold arcs. These occur in pairs with a common tangent at the origin, emanating in the same 
or opposite directions from the origin according as m is even or odd respectively. In the even case, 
as e crosses the first arc in the appropriate direction a pair of solutions is created away from S 
by a saddle-node (fold) bifurcation, and as e crosses the second arc a second saddle-node pair of 
solutions is created close to the n-images of the first pair. In the odd case a saddle-node bifurcation 
occurs across each arc. 

End arcs (m even). These occur in pairs with a common tangent at the origin, emanating in 
opposite directions from the origin. As e crosses an end arc in the appropriate direction a pair of 
solutions is created away from S in a saddle-node bifurcation; the two directions are in opposing 
senses of rotation about the origin. 

Hysteresis arcs (m odd). These occur in pairs with a common tangent at the origin, emanating 
in the same direction from the origin. As e crosses the first arc in the appropriate direction a saddle- 
node bifurcation creates a pair of solutions, one of which coalesces with an existing solution (via a 
hysteresis loop) in a saddle-node bifurcation across the second arc. 

This elaborates on and gives a geometric setting for Theorem 2. IB of [23]. See Figure [TTj and 
compare with Figures 3 and 4 of Note that the interpretation regarding solutions {x,y) 

passing through y — Q does not hold in the general case m > 2 since clq — Q does not imply — Q 
for points on A^f . 



Remark. These results on the structure of A/^ and the form of the expression (4.9 1 show that for 
£ g E and d—k=l the generic solution branches are of the form {e,x,y) ~ (i™,i™+^,i). This 
applies also in the case m — I (normal nondegeneracy) where y represents a coordinate transverse 
to S and the nonzero solution value y varies linearly with e. 

4-. 3. The variational case 



As noted in Section 2.1 when the original vector field is the gradient of a smooth real- valued 
function then the reduced vector field may be regarded as the gradient of a real- valued function 
/g . Therefore in this setting the bifurcations of equilibria are controlled not as in Section [s] by IC- 
versal deformation of the map germ 

Fq{xq, •) : M*^ M'* X R'' 

at e M'' but rather by an 7?.-versal deformation of the function germ 

/o(a;o, •) : M'^ ^ M 

at G M*^ where Fq = grad/p. In the case of uniform normal degeneracy of S the germ fo is 
independent of xq e S. 



To study the bifurcation geometry we follow the formal procedure of Section 3.2 but now 
in addition to the simplification in dealing with functions rather than vector fields there is an 
added twist: the gradient of a deformation of fo{xo, •) will involve derivatives with respect to the 
'parameter' x as well as with respect to the variable y. 
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Figure 11. Schematic bifurcation diagram for d = k = l,q = 2 and m odd (see text), showing 
hysteresis arcs (unbroken lines) and fold arcs (dashed lines) through the origin in parameter space 
R^. Numbers indicate the number of solutions in each region. Note that for each opposite pair 
of fold arcs the two folds approach from opposite sides in the plane. 



Assume we are in the uniform case, and write fo{x,y) = fa{y). If /o has an r-parameter 
7?,-versal deformation of the form 

aivi{y) H h arVr{y) + fo{y) 



then as in Corollary 3.9 we express as 

fe = ai{e, x)vi{y) H h ar{e, x)vr{y) + fo{y) 

and now (after dropping the tilde) we find = {dxfe, dyf^) where 
dxfe = a[{e,x)vi{y) H h a'^{e, x)vr{y) 

dyfe = ai{e,x)v[{y) H h ar{e,x)v'^{y) + Fo{y) 

and we write a[,v'j for grad^Oi, gra,dyVj respectively. The condition for {s,x) to belong to 
the discriminant is an algebraic condition on the 1-jet (with respect to x) of the map 



(4.18) 
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a = (ai, . . . ,ar) : M'^ x S* ~> W . We illustrate this with a simple but non-trivial example that 
exhibits both the similarities with and the differences from the general geometric formalism of 
previous sections. 

Example. Let fc = 1 and consider the bifurcation problem ~Q where = grad with 

An 7?.-versal deformation of the function germ is 

ao + aiy + • • • + + y™ 

(see [6],|25]) and so for to = 3 we can write 

fe{x, y) = ao(£, x) + ai(e, x)y + y'^ 



and we have 



dxfe{x, y) = ao(£, a;) + a[{e, x)y 
dyfe{x, y) = ai{e, x) + 3?/^. 



The condition for (e,x) to belong to the discriminant Aj? in the case d = 1 (compare (4.3 1) is then 

3(a^(e,x))2 + {a\{£,x)Yai{e,x) = 0. 
Following the geometric analysis of Section [4] we write e — ps in polar coordinates and 

a{ps,x) ^ ph{s,x) + 0{p^) (4.19) 
and find that the conditions for &(e, x) to lie in the positive tangent cone T are 

6o(a;) = 0, hi{s,x)<Q. (4.20) 

where 6q = grad^ ho- For simplicity of description we take as usual the case g = 2 so (s, a;) £ x S". 
Generically the locus 

B'a = {(s,x) G §^ X 5 : 6o(s>2:) = 0} 

is a collection of smooth curves (codimension-1 manifold) such that the singularities of the projection 
p\B'i^ : B'q ^ §1^ are isolated folds. There are at least two of these on each compact connected 
component of S . (For g < 4 the singularities of p will generically be folds, cusps {q — 3,4) or 
swallowtails {q — 4), while for g > 4 singularities of corank > 1 can appear. See [3BJ Ch.VII,§6] for 
example.) Likewise the locus 

Bi = {(s,x) e §^ X 5 : bi{s,x) = 0} 
is generically a smooth curve (codimension-1 manifold) which intersects Bq transversely at points 



which are not singular points of the projection p\Bq. The analysis proceeds as in Sections 4.1.2 



and 4.1.3 noting that now end-points are given by Bq n Bi while intersection points are the critical 



points of 6i|i?o • Thus the generic bifurcation geometry turns out to be the same as in Theorem 4.3 
for the (non- variational) case to = 2 although the interpretation of key features is here tied to the 
variational structure. 

For the general case fo{y) = we define the 1-jet map 

/a ■.R'^ X S ^ J^(S',M'"~^) : {e,x) (x; ao, ai, 02, . . . , a,„_2; Oq, a'l, . . • ,a'„_2) 
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where — ai{s,x) and — grad^ 
consider the map 



for i — 0,1, 



2. Taking local coordinates on S we 



1) - 1 (see 

^-i) ^ E-* X K 

m— 1 



y^(0,...,0,my™-i) 

which has /C-codimension e := (d + l)(m - 
with discriminant Ah C M"^. Locally J"'^(S'. 
let Ay = (tt o jia)^^AH where tt is projection of J^{S. 
omitting x and ao). We then consider the related map 

b : §9-1 X S -^R" : {s,x)^ {bi,b 
-1 X S* ^ I 



) and a /C-miniversal deformation 

,m-l ^ ]gd(m-l) = X ]ge+l 

) onto the last e coordinates (that is, 
,b'o,b[,...,b'^_,) 



with 5 : §''-1 X S ^ M'""^ defined as in (4.191. Under suitable non-degeneracy hypotheses on b 
which, using the Thom jet transversality theorem (see [6], [22], [25], [26]), we expect to be generically 
satisfied by the map , local geometry of the discriminant and its projection to the parameter 
space M'? (giving the bifurcation structure) can in principle be determined — at least for small 
values of q and m. We leave details to the enthusiastic reader. 



4-.4- Possible further reduction 

In the general case a further reduction is technically possible, using the IFT one more time. For 
simplicity of illustration take g = 1 so that e G M. 

Given x £ S, write g{x) — T{0, x, 0) where T : RxU ^ M'^+'^ is the deformation map as in (2.51 
with (x, y) G U . We make the following assumption, which holds generically as dimS" = d < d + k: 



(A4) The image oi g : S —> R''-^^ avoids the origin in 

This means that given xq £ S there is at least one component gi of g that does not vanish at (xq, 0) 
when e = 0. By the IFT the component of = can then be solved close to u = (xq, 0) in 
the form e — £i{x, y) . Substituting this into the other components of yields d + k — 1 equations 
to be solved for {x, y) close to (0, 0) G M"* x M*^. The solution set can be viewed as the zero set for 
a d-parameter family of smooth map germs R^ E''+'^~i, where techniques of versal deformation 
still apply, now in one lower dimension in parameter space and target. 

For example, if d = 1 then we have a 1-parameter deformation of germs R^ M*^ about which 
much is known for small k: see [35j for the case k — 2. li d — 2 then we are dealing with germs 
_^ ^k+i. ^ fQj. i^jjg (,g^gg ^ ^ 2 here. 

In this formal reduction of the problem, however, the controlling role of e has been sacrificed, 
and it is not clear that there is practical benefit to understanding the bifurcation structure. We do 
not pursue this approach further. 

4-. 5. Nontrivial bundles 

In the case of nontrivial bundles N S,TS,K etc, where the same coordinates {x,y) cannot be used 
on the whole oiU <Z K, the deformation theory must be re-expressed in terms of jet bundles. Given 
smooth vector bundles D,E over S and a family of bundle maps : D ^ E with e G M'', the 
generic structure of F^ is expressed in terms of transversality of the r-jet extension 

fF -.R"^ xD^ J''{D,E) : {e,u) ^ fF,{u) 
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to suitable stratifications of the r-jet bundle J^{D,E) and for sufficiently large r (cf. [5 ). In our 
context this theory has to be applied to the reduced vector field : K ^ TS P. The geometry 
of the zero set is captured in local coordinates using /C-versal deformation theory as we have done, 
but global aspects are dictated by the geometry of the bundles. It would be interesting to explore 
this further even in the simplest case c? = fc = 1 with K both the trivial and nontrivial M-bundle 
over a circle. 



5. Local branching analysis in the case q = k = I. 

It is one thing to understand generic structure of the bifurcation set, but another to recognise the 
specific branching geometry in particular cases. In this section we focus on the latter. We return 
to the key case q = k — 1 , taking advantage of the fact that a singular function of one variable 
2/ S K is dominated by the first nonvanishing term of its Taylor series in order to classify the local 
branching behaviour of Fg in a systematic algebraic way (bottom up) that complements the generic 
geometrical approach (top down) that we have described so far. 



Recall from (2.51 that we have in local coordinates 

F{e,u) = Fo{u) + e^{e,u) (5.1) 

where u — {x,y) G S x M. and e e M. Each component of Fq{u) has a Taylor series in y 
with coefficients depending smoothly on x G S. We do not immediately assume here that the 
normal degeneracy is uniform, but we suppose some control on degeneracy by making the following 
assumption: 

(A5) Each component of Fq(u) has nonzero Taylor series in y at every point x G S. 

In other words, the normal degeneracy in the problem is of finite order everywhere on S . This 
is not essential in all that follows, for it will be clear how to amend the discussion to allow infinite 
degeneracy in one component. We leave aside, however, the complications of infinite degeneracy in 
more than one component. 

5.1. Necessary conditions for branching 



From (2.3 1 and (2.4 1 the i^^ component of F{e,u) can be written in the form 

F,{e,u) = eg^{e,u)+y'^'r,{u) (5.2) 

for some integer > 2 and a smooth function : t/ — s- K ; we suppose that ri(-, 0) is not identically 
zero, otherwise we would increase the exponent . Let Xq G S and number the coordinates so 
that at xq 

mi < m2 < ■ ■ ■ < rud+i- 

Proposition 5.1 Necessary conditions for xq G S to be a branch point are that for 1 < i < J < d+1 
fidj = when mi < mj 

(5.3) 

dij = when mi = mj 
for aUl<i<j<d+l, where gi — gi{0, Xq, 0) , rj — rj(xo, 0) and dij :— rtgj — rjgi . 
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Proof. From ( |5.2[ ) we see that 

F,{e,u)gj{e,u) - Fj{e,u)g^{e,u) = y"''r^{u)gj{e,u) ~ y"'^rj{u)g^{e,u). (5.4) 
Dividing through by y™' and letting (e, x, y) (0, xq, 0) immediately gives the results. □ 



The conditions (5.3 1 are not all independent. Given that the assumption (A5) holds, choose 
xq E S and let 

n — n(xo) — max{i : g,; ^ 0} . 



Then the conditions (5.3) reduce to 



(a) Tj = when > m„ 

(5.5) 

{h) dij = when nii =^ rrij = m„ . 



Let Jn = {j '. rnj = m„}: then J„ is an integer interval say, with I < n < p. In (5.5 1 the 
number of conditions (a) is Z — 1 and the number of independent conditions (b) is p — I. Thus (5.5 1 
represents p — 1 conditions, which together with 

9n+i = ■ ■ ■ = gd+i = 
give a total of d — n + p conditions on xq G S. Hence we have the following result. 



Corollary 5.2 Generically no points xq £ S satisfy the conditions (5.5) when n < p, and there is 
a set of isolated points xq € S which satisfy (5.5) when n — p. 

In the uniform case the results can be expressed more simply. Each r,j is nonzero everywhere on S 
and the exponents do not depend on x G S' , so if branching is to occur at any xq S S* we must 
have n — p (since if n < i < p then d„i — implies = 0) and also I — 1 (else rj = for j < I). 
Hence we obtain in this context the neater result: 

Corollary 5.3 In the uniform case a necessary condition for Xq ^ S to he a branch point is that 
9[n]{xo) = where 

9[n\ ■= {d2, ■ ■ ■ , dn,gn+i, . . . , ffrf+l) : S ^ R"^ 
in which dj denotes dij for j — n, . . . , d and n = max{i : = mi}. 

Of course, if nii < m2 (so mi < nij for all j > I) then g[„] contains no dj terms and the result is 
even more straightforward : 

Corollary 5.4 // mi < m2 then branch points xq S must be zeros of the map 

g:^ {g2,...,gd+i):S^R''. 

Generically € M'' will be a regular value of g and the zeros of g will be a discrete subset of S. In 
any given context there may be analytical or topological tools (for example, index theory) available 
to locate or count them. 



The previous results gave necessary conditions for solution branches to emanate from 
(0, xq, 0) e M X [/ c M X 5 X M. We now turn to sufficient conditions. 
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5.2. Sufficient conditions for branching 



Let xq d S satisfy the necessary conditions (5.5) for branching. With the notation of Sect ion 5.1 



since (?„ 7^ at (0,a;o,0) we can use the IFT to solve the equation Fn{e,u) = from (5.2 1 in the 
form 

£ = £^{u) = y'^"S,,{u) (5.6) 

where (S„ is a smooth function on a neighbourhood of (a;o, 0) in C/ C K. Substituting this into each 
component Fi with i =/= n and dividing through by appropriate powers of y we see that finding a 
solution branch amounts to finding solutions u — (x, y) close to (a;o, 0) to the set of equations: 

Sn{u)gi{en{u),u) + y"''-"'"r,{u) = 0, i>n 
dnj{en{u),u) ^ 0, l<j<n 

y"'-'"'^S^{u)g,{£„{u),u) + r,{u) = 0, j<l. 

Again by the IFT, these equations will have a smooth solution branch in the form x = x{y) (that 
is, transverse to Sq in U) if and only if the Jacobian matrix of the map 

S -.x^ (ri,...,r;_i,d;,...,d„_i,(5„5„+i,...,(5„5d+i) (5.7) 

(there are no dj terms if n = ^ ) is nonsingular at Xq . Now suppose 

r„ :=r„(xo,0) ^0 (5.8) 

which is satisfied automatically in the uniform case. Then solving Fn{e,u) = shows that 
Sn{xn,0) ^ 0. Difi^erentiating Sngi at {0,xo,0) gives 

da;{Sngi)iO,xo,0) = S„{xo,0)dxgi{0,XQ,0) 

for z > n as gi{xQ,0) — 0. Therefore, since Sn{xQ,0) =/= 0, the Jacobian condition above can be 
re-stated as follows. 



Proposition 5.5 Given that the necessary conditions (5.3) are satisfied at xq Cz S and (5.8) also 
holds, a sufficient condition for the existence of a solution branch emanating from (xqjO) and 
transverse to So in K is that xq be a regular zero of the map 

gi,n := {ri, . ■ . ,r/„i, d;, . . . ,dn-i,gn+i, ■ ■ ■ ,gd+i) : S —^W^, 
the r and d terms being absent when n = I. 

In the uniform case (where all rj =/= 0) the expression is again much simpler: since mi = . . • = m„ 
and there are no rj terms we have: 

Corollary 5.6 In the uniform case, a sufficient condition for xq (z S to be a branch point is that 
Xq be a regular zero of the map g[n] ■ S — > M''. 

As before, note that if mi < m2 then g^^^ becomes the map g : {g2, ■ . ■ ,gd+i) ■ S — > M.^. 

Remark. We can regard the normally nondegenerate case as an extension of the above: here 
mi = 1 and all other m^ > 1, and so regular zeros of g are branch points. 

Remark. In general if a zero xq of gi^n is not a regular zero then progress can still be made by 
standard bifurcation techniques. For example, if the Jacobian matrix at xq has rank d — 1 then a 
further Lyapunov-Schmidt reduction leads to a single bifurcation function of one variable, whose 
zero set can be determined locally from the quadratic terms if nondegenerate or from higher order 
terms by standard singularity-theory methods otherwise. 
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5.3. Examples 



We illustrate the preceding formal descriptions by some simple examples. 



Example 1 (dim 5 = 1) 




From Proposition |5.1| a necessary condition for G 5 to be a branch point is that xq be a 
zero of the function 



X i-> ri(a;, 0)32(0, x,0) 
and a sufficient condition is that xq be a simple zero. 



Example 2 (dimS" = 1) 



Fa{x,y) = {ri{x,y)y'^,r2{x,y)y'^). 



With (51,52) as in (5.9 1 a necessary condition now is that xq be a zero of the function 



di2 : X 51(0, X, 0)r2ix, 0) - 52(0, x, 0)ri{x, 0) 



and a sufficient condition is that xq be a simple zero. 

In the uniform case the functions ri, r2 are nonzero constants and (as we have seen in general) 
the necessary and sufficient conditions are conditions on (71,(72 alone. Generically 52 (Example 1) 
or di2 (Example 2) will have a discrete set of zeros on S , all of them simple. 



Here with the usual notation we see that necessary conditions for G 5 to be a branch point 
are that rig2 — rig^ = r2g3 = at (x*, 0) G x M, or in other words that a;* is a zero either 
of (^1(72, 33) or of (ri, r2), these maps S M.'^ being evaluated with y = 0. In the uniform case 
this reduces to being a zero of (52,53); then being a regular zero is a sufficient condition for 
x^, to be a branch point. 

We leave to the reader the extension of this example to the cases when some components of 
Fq have leading terms with equal powers of y . 

5.4. The variational case 

We now consider how this local branching analysis applies in the variational case. For simplicity 
we take q = k = 1. 

Following the preceding sections and dropping the tilde notation let / : R x [/ — s- R be a smooth 
function such that dfo{u) = and d^fo{u) — for all u = (a;,0) G So, where /^(u) = f{e,u). 
Assuming the non-flatness condition (A5) we write / in the form 



Example 3 (dim 5 



FQ{x,y) = {ri{x,y)y'^,r2{x,y)y^,r3{x,y)y'^). 



/(e, X, y) = £5(e, x, y) + y"''r{x, y) 
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with m > 3 and r(-,0) not identically zero on S; in the uniform case r{-,0) ^ on 5. Critical 
points of /e : [/ ^ M occur where 

r e d^g{e, x, y) + 9,r(x, y) = G M'' 



y e dyg{e,x,y) + my"^ ^r{x,y) + y^"- dyr{x,y) = £R. 
Proposition 5.7 A necessary condition for Xq Cz S to be a branch point is 
r{xo,0)d.,g{0,xo,0) = OeR''. 

leads (after cancelling y™~^) to the equation 



Proof. The argument of Proposition 



5.1 



m r{x, y)d^g{e, x,y) + y [d^g{e, x, y) dyr{x, y) - dyg{e, x, y) d^r{x, y)] = 



and as (e, x, y) — > (0, xq, 0) we obtain (5.11 1 



Corollary 5.8 In the uniform case all branch points are critical points of g(0, -,0) : Sq 



(5.10) 



(5.11) 



□ 
□ 



Proposition 5.9 If {5.11) holds for xq € S, then a sufficient condition for xq to be a branch point 
is provided by either of the following conditions at xq : 

(i) dr/dxn^i 7^ and xq is a regular zero of the map 

X ^ {mr, 62, ... , e„_i, /i„, . . . ,hd) : S ^ 



where e,- 



-^.^ and h, = 4^: 



'^3 ' dy dxj- _ 

(ii) r 7^ and xo is a nondegenerate critical point of g{0, - ,0) 



Proof. If dj:g{0,xo, 0)^0 then r{xo, 0) = by ( 5.11[ ). Condition (5.8) here is dr/dxn-i ^ 0, where 



(with the notation of Section 5.1 1 we have J„ = [2,d+ 1] and (i) follows as in Proposition 5.5 
that dj reduces to ej since r 



0. 



note 

□ 
□ 



If dxg{0, xa,0) = then again from Proposition 5.5 with (5.8 1 we obtain (ii). 

Corollary 5.10 In the uniform case r ^ and it is only the condition (ii) that applies. 

Here, just as in the normally nondegenerate case, we can use topological tools such as Lyusternik- 
Schnirelmann category or Morse Theory in order to give a lower bounds for the number of branch 
points: see [3j[9l|35]. Observe that this global result characterising the branch points is the same 
as in the normally nondegenerate case. However, the branching behaviour itself is different: recall 
the Remark following Theorem |4.6[ 



6. Applications 

Here we outline some examples of problems where the techniques presented in this paper can be 
applied. 
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6.1. Periodic orbits 

A standard technique for finding periodic orbits (of period T) for a vector field Xq is to consider 
them as zeros of the nonlinear operator 



defined on a suitable Banach space Ti. of T-periodic functions. There is a natural action of the circle 
group on H given by time-translation modulo T, and if Fq{x) — then also Fq{xs) — for all 
s G = M mod T, where Xs{t) denotes x{s + t). Thus each zero of Fq is automatically a member of 
an S-'^-orbit S of zeros. When this orbit is normally nondegenerate the problem essentially reduces 
to a problem on S and techniques that are now standard (jH [131 [HI [Ml 1211 132]) can be applied 
to study bifurcations of zeros of Fq under perturbations of Xq . The first authors apparently to 
attempt a systematic approach to normal degeneracy in this context were Hale and Taboas [IS] 
who used analytic methods to study a differential equation of the form 

X + g{x) = eih{t)x + e2f{t) e M 

close to a degenerate T-periodic solution of a; -f g{x) — 0, where h and / are T-periodic functions. 
As indicated in Section |4.1.3| our Theorem |4.3| recovers their results and sets them in a geometric 
context. 



6.2. Hamiltonian systems 

In the case of Hamiltonian systems, where periodic orbits may typically lie in flow-invariant tori, 
each zero of Fq will lie in a larger-dimensional manifold (torus) of zeros. We expect normal non- 
degeneracy of this manifold within each energy level, but for discrete values of the Hamiltonian 
function (energy) we may find normal degeneracy. This appears both in the study of the Poincare 
map for the periodic orbits and equivalently in the study of manifolds of critical points of the Jacobi 
metric. 

To illustrate this in a relatively simple situation where the periodic orbits are isolated in each 
energy level, consider a Hamiltonian system in M? with a smooth Hamiltonian of the form 

H{q,p) = K{p) + V{q) 

where (q,p) G T*M? ^ x , and where K{j)) — |pp/2. Up to parametrisation, a periodic orbit 
q{-) for the system is a critical point of the Jacobi metric (see |T] for example) 

J[q]^ f {E-V{q{s)))K{q{s))ds 
Jo 

which has natural time-translation § -invariance as a smooth function on H = i?^([0, 1],M2); here 
E is a. constant greater than V{q) throughout the region of under consideration. 

Assume that the system is rotationally symmetric about the origin, so that V{q) — V{r) with 
r = The Euler-Lagrange equation derived from is 

±[q{E-V{r)]+^-^^-V'{r)q{t) = (6.1) 
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where ' denotes J^. Now let qo{t) be a uniform circular orbit written as qo{t) = e(i) where 
e{t) e is a unit vector. Substituting <7o(^) into (6.11 and using the identity e(t) = — |e(i)pe(t) 
gives 

{E~Viro))^'^V'iro) (6.2) 

as a necessary and sufficient condition for qo{t) to be a periodic orbit for the system. Note that 
circular periodic orbits of fixed radius form a manifold Se in Ti. that is a copy of S^: 

SE = {qseH:qs{t) = Rsqo{t)} 

where Rs G 5*0(2) denotes rotation through angle s G [0,2tt). 

To ascertain normal (non)degeneracy we next study the second derivative D'^J'[qo] which is 
defined on on the space H and corresponds to a den sely defined self-adjoint operator on i^([0, 1], M^). 
We show in the Appendix (Section 7. 1 ) that given ( 6.2 1 the kernel of D^J'[qo] has dimension 2 when 

roV"{ro) = {2n-' -9)r{ro) (6.3) 



for some n € Z, and has dimension 1 (corresponding to Tq^S) otherwise. Note that condition (6.3 1 
can be read in terms of discrete values {En} of the energy E : 

^" = ^^ + 2(2;^^"^. (6.4) 

For energy values En the corresponding Se„ is a normally degenerate manifold of points representing 
periodic orbits. 

A potential V{q) for which all this can be explicitly verified is the so-called Mexican hat: 
V{q)^-^\q\' + \\q\'. 

The above methods allow us to detect the presence of periodic orbits that persist after perturbing 
V in the annular Hill's region given by {g e M2 . yj-^-) < ^1 ^Yiem - ^X^ < E <0. Setting E = En 
and fixing u G Kq^ where Kq^ is a complement to Tq^SE,^ in k.er D^J^lqo] it is not difficult to check 
that the Taylor expansion about qo of the Jacobi metric after Lyapunov-Schmidt reduction reads 

Jiy) J[qo + yu]= J[qo] + Joy' + 0{y') (6.5) 

where Jo = jfD'^J'[qo]{u,u,u,u). This places us in the context of Section 
m = 4. 

An example of bifurcation can be seen by considering a potential 
V,{q) = V{q) + EiMq) + £2M<l) 

where (t>i{q) — (j)i{qi,q2) is a smooth function, i — 1,2. Up to degree 4 in y the reduced Jacobi 
metric is now 

j,[qo + yM] - JM + Joy'* + ei«'i[go,y] + £2«'2[go,y] (6.6) 

where 

*»[90,y] = - / dt\\q^{t)f{(^,{qQ{t) + yu{i)) 
Jo 

for i = 1, 2. Now let qo{t) — (tq cos(ct; t + x),rQ sin(ti; t + x)) be an element of Se„ parametrised by 
X ^ S^. The expression ( 6.6l reduces to 

J{e,x,y) J[go] +Joy'^ + £igi{x,y) + e2g2{x,y) 



4.3 



with d = k = I and 
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up to degree 4 in y, where 

gi{x,y) = '^i[qo,y] 

for i = 1,2. Following the Example given in Section |4.3| but now with m = 4 we see that the 
bifurcation geometry is determined by the geometry of the curves B'q and Bi in §^ x S* given by 
6q(s,x) = and bi{s,x) = respectively, where 6o;^i £^re derived from the perturbation g as in 
Section l4?2l with m — 4. 



Remark. The case of a single parameter s (so q = I) can be interpreted in this context as a fixed 



choice of s € Here the function g(0, -,0) : 5*0 ^ M of Section 5.4 corresponds to the function 



5o(s, •) of Sections 4.3 and 4.2 The necessary condition for xq to be a branch point as given by 



Proposition 5.7 corresponds to the statement that the curve B'q intersects the circle {s} x 5 at 



{s,xo), and the sufficient condition given by Proposition 5.9 li) corresponds to the statement that 
it does so transversely. Since every smooth function on a circle has at least two critical points there 
will be at least two such intersections (Lyusternik-Schnirclmann category) , and in any case an even 
number if they are transverse (Morse theory). 

6.3. Steady states in chemical reaction networks 

In kinetic models for chemical reaction networks the time evolution of a vector x(t) — 
{xi{t), ...,Xn{t)) of concentrations of n chemical species involved in r reactions is given by 

dx 

— = F{x,k) -.^ Biy{x,k) (6.7) 
where i? is an n x r matrix (stochiometric coefficients) and 



is a smooth map, often in fact a polynomial map. A review of this type of system can be found 
in p3 . The components of fc g M"^ are the reaction constants. For given k the stationary states of 



(6.7 1 are thus given by the set 

Sk := {x e M" : Vk{x) £ ker(S)} 
where Vkix) = ^{x, k). 

In applications, one is interested in studying steady states of perturbed systems of the form 

F^{x) B v{x, k) + e0(£, x) (6.8) 

for e e M"? where 0(e, x) is smooth (not necessarily polynomial) function. It is therefore appropriate 
to apply our analysis to systems of this type. 

Let R C M" denote the range of B and let tt : K" ^ i? be a projection onto R. Write 
B = -K o B : W ^ R. If G i? is a regular value oi Bov^. (that is, the map Vk is transverse to ker B) 
then Sk is a smooth submanifold of K" with codimension equal to the dimension of R (the rank of 
_B); moreover, dimS'^ = dimker £'i^(2:) for x ^ Sk &o the manifold Sk is normally nondegenerate. 

If is not a regular value oi B ovk it means that there is at least one point a; € S'fc at which Vk 
fails to be transverse to ker B. Situations can arise where this occurs simultaneously at all points 
of Sk and we have uniform normal degeneracy. We present now a simple example of this class of 
system. 
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Consider a set of reactions of the form 

1 ,^1 



Xi + X2 — > , X^ — > Xi + , — > X2 

where the first two reaction rates are fixed to 1 and the third is a function v of the concentration 
[Xi] of Xi. The third reaction is interpreted as an external input whose rate is a function of 
[Xi]. Now let [Xi] = Xi denote the i-th concentration, so a; = {xi,X2,X3) G (all components 
non-negative) and the dynamic equations read 

Xi{t) = -Xi{t) X2{t) + X3{t) 



X2{t) = -Xi{t) X2{t) + X3{t) + V{xi{t)) 



(6.9) 



is given by 



isit) = xi{t) X2{t) - X3{t) 
where here we take k as fixed and omit it from the notation. The map ly : M.^ 

v{x) = {xiX2 , v{xi) ,xz)^ , 
and the matrix B is 

/ -1 1 
B= -1 1 1 
\ 1 -1 ^ 

with 1-dimensional kernel ker B = span{(l, 0, 1)"^} and range R = span{(0, 1, 0)"^ 

Assume that the rate v(-) is a smooth function with a unique zero xl > 0. 
states are therefore 

5*+ = {{xi,X2, X3) : X1X2 = A, v{xi) = 0, 2:3 = A} 
for A > 0, which is the part of the 1-dimensional afline space 

S = {{xi,X2,X3) : xlx2 - X3 ^ 0, xi^xl} 
where all components are positive. The tangent space of S* at x S S' is independent of x, namely 

T,5 = span{(0,l,a;I)^}. 



(1,1, -in- 

The stationary 



We have 



and 



Dv{x) 



BDv{x) 




-1 



Hence we see 

Proposition 6.1 The map Bv :M? ^ R has as a regular value if and only ifv'{x*i) 7^ 0. □ 

Now suppose that v'{x\) — 0. This does not alter 5* but does destroy the transversality of v to 
ker B. If 2; S 5*+ C S then T^S is annihilated by BDv^x) but not by Dv{x), so that 5+ is normally 
degenerate with constant corank 1. We choose a suitable complement to T^S in ker i3I?zy(a;). 
In the Appendix Section 



7.2 



it is shown that the reduced equation Fq{u) — (see Section |2| has 
uniform quadratic degeneracy provided v"{x\) ^ 0. For studying bifurcations we are therefore in 
the context of Section WA\ with d = k = 1 and m = 2. 
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7. Appendix: calculations for examples in Section [6} 

7.1. Hamiltonian periodic orbits 



In this section we derive the condition (6.3). 

Let q{t) — r e{t) G be a circular path with constant speed \e{t)\ = ui. Let u{t),v{t) G H. 
Direct calculation shows that the Hessian oi J a± q is 

D^J[q]{u,v)= [ dt[{u{t),i{t)){E-V{r))] + 
Jo 

-uo f dt [{v{t),e{t)){u{t),e^{t)) + {u{t) , e{t)) m , {t))] rV {r)+ 



to 



2 



dt 







^ ( ^"W + -V'{r) ] {e{t)Mt)) {e{t),v{t))+rV'{r){u{t)Mt)) 



To compute D J aX q — q^ we set r — rt^ and use condition (6.2 1. We find 

D^J[qa]{u,v)^ I dt[{u{t),v{t)){E-V{r^))] + 



-2c. / dt [{{v{t),e{t)){ii{t),e^{t)) + {u{t),e{t)) m,e^{t))) {E - V{t^))] 
'I 



uj'' I dt 



{rlV"{r,) + 2{E- V{ro))) {e{t),u{t)) {e{t)^v{t)) + {E - V{ro)) {v{t),u{t)) 



Since E ~ V{ro) > the condition D'^J'[qo]{u,v) = is equivalent to 

dt [{v{t),u{t))^2uj{{v{t),e{t)){u{t),e^{t)) + (u{t),e{t)){v{t),e^{t))) + 

n {e{t),u{t)) {e{t),v{t)) - {v{t),u{t))] = 



where 



1 



(7.1) 



(7.2) 



^2(£;- V^(ro)) 

Integration by parts and a standard argument from calculus of variations shows that u{t) lies 
in the kernel of D^j7[(?o] precisely when 



u{t) + 2io e{t) {u{t),e^{t)) - 2lo -{{u{t),e{t))) e^{t)- 



+2^2 {u{t),e{t)) e{t) + Y {e{t),u{t)) e{t) + u{t) = 0. 



(7.3) 



Decomposing u{t) as 

u{t) = a{t) e{t) + b{t) e^{t) 

we find the equation (7.3 1 splits into 

'd{t) + {A + n/2)uj^ a{t) = 0, 
bit) = 0. 



(7.4) 
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Equations (7.4 1 imply that 

dimkcT D^J[qo] = 



1 if (4 + n/2) ^ with n e Z, 

2 if (4 + n/2) = for some n G Z. 



(7.5) 



Finally, the expression (7.2) for H. and condition (6.2 1 show that the condition (4 + n/2) = n? is 
equivalent to (6.3 1. □ 

7.2. The reduced equation for the chemical reaction network 



Here we construct the reduced equation for (6.9 1 perturbed as in (6.8 1. From (6.9 1 we write (6.8 1 
as 



F^{x) := Bv{x) + x) = f{x) w + v{xi) e + e(f>{e, x) 



(7.6) 



where w — (1,1,-1)^ 



(0, 1,0)"^ and f{x) = —xiX2 + 2:3 £ M. For x £ S+ the range of 



(7.7) 



e 

BDiy{x) is spanjii;, e} except if w'(xj) = in which case — spanjw}. For z G we have 
{w,F,{z)) = f{z)\\w\\^ + v{zi){w,e) + {w,(p{e,z)) 

= 3/(z) + v{zi) + {w, (j){e, z)). 

For x G Sj\- and v' {x\) = we have ker BDv{x) = where u ~ {x2^ a^i, —1)"^ and so it is natural 
for the Lyapunov-Schmidt reduction to take Lx — span{u}. Accordingly for 2 in a neighbourhood 
of S+ in we write 

z = a; + j/n + Au G © 

where x G 5"+ and 

Kx — {T^S ® span{u}) = span{n} 

where n = (1 + x*i^ , —sx\, s)^. By the IFT we can solve 

{w,F^{x + yn+ \u)) (7.8) 

in the form A = X*{x,y,e) where A* is a smooth function defined on a neighbourhood of S'o x in 
X M'J and with 7^(2^, 0, 0) = for all x £ S. It then remains to solve the reduced equation pair 

{s,F^{x + yn + X*{x,y,e)u)) = 0, {p,Fs{x + yn + X*{x,y,e)u)) = (7.9) 

where spanjs} — T^S and span{p} = Pj. chosen so that does not lie in R^ ® T^S (for example, 
p= (0,1,0)^). Now 

{s,F^{z)) = f{z){s,w) + v{zi){s,e) +e(s,0(e,z)) 



where 



a(s) 



= a{s)v{zi) + £(6(s), 0(6, z)) 
1 , 



[S,W) 



and 



1 



6(s) = s — -(s, w)w 
o 



(7.10) 



from (7.7l, and similarly for {p, Fi^(z)). With z = x + yn + X*u we can write 

v{zi) = v{xi + yni + X*ui) — v{xl + y{l + xl^) + X*X2) 

= i«"(xt)(y(l + xf) + X*uif + 0(3) 

where 0(3) terms are of third or higher order in {yi,e) and smoothly parametrised by x G S+. 
Since A* has no linear term in y the smooth map yi y{l + Xi'^) + X*ui is a diffeomorphism germ 
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at the origin in = M smoothly parametrised hy x E S+ and sufficiently small e e M''. After 
we make this change of coordinate, the assumption v"{x\) ^ allows us to apply and a further 
local diffeomorphism in to absorb the 0(3) terms. We may therefore suppose without loss of 
generality that the equations ( 7.9 1 are 

a{s)cyl 
a{p)cyl 



- e{b{s),(l){e,x,y)) 
' e{b{p),(t){e,x,y)) 



6.3 



More 



where c = ^v"{xl) ^ 0, and so we have a quadratic system as claimed in Section 
specifically, the problem is of the form given in Example 2 in Section |5.3| and hence the branch 
points X £ S+ are given by the zeros (provided they are simple) of the function 

g{x) := a(s)(&(p),0(O,x,O)) -6(s)(a(p),0(O,x,O)). 
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